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Invited  Presentation 


Army  Tank-Automotive  Propulsion  Activities 


W.  Bryzik 

US  Army  Tank-Automotive  RDE  Center  (TARDEC) 
Warren,  MI  48397-5000 


The  propulsion  system  forms  a  critical  building  block  for  present  and  future  military  ground 
vehicles.  Characteristics  of  particular  importance  to  TARDEC's  ground  propulsion  systems 
include:  low  system  volume  (including  engine,  cooling  system,  air  handling  system, 
turbomachinery,  ducts,  etc.),  low  system  weight,  improved  fuel  economy  or  vehicle  range,  rapid 
acceleration,  reduced  signature  (including  thermal  acoustic,  visual/smoke,  etc.),  high  reliability  and 
survivability,  and  improved  cost  effectiveness.  Innovative  integration  of  the  propulsion  system 
components,  when  combined  with  individual  component  technology  gains,  produce  an  exceptional 
opportunity  for  optimizing  the  above  listed  characteristics. 

Various  propulsion  system  activities  at  TARDEC  are  discussed  during  this  briefing  from  an 
overall  perspective  standpoint.  Representative  example  projects  dealing  with:  1)  low  heat  rejection 
combustion,  2)  cold  start  characteristics,  3)  engine  finite  element  analysis  procedures,  4)  advanced 
diesel  technology,  5)  engine  wear  analysis  and  life  prediction,  6)  thin  ceramic  coating  performance 
optimization,  7)  multi-fuel  adaptation  to  two-stroke  engines,  and  8)  a  number  of  other  efforts 
aimed  at  optimizing  the  characteristics  given  above  are  discussed.  In  addition,  the  importance  of 
TARDEC's  leveraging  of  propulsion  activities  being  funded  by  sources  such  as  the  Army  Research 
Office  (ARO),  Army  Research  Laboratory  (ARL),  ARPA,  DOE,  NSF,  and  others  is  noted  with  a 
few  examples  given. 
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Invited  Presentation 


Vehicle  Propulsion  Directorate  Engine, 
Combustion  and  Reacting  Flow  Interests 


R.  Bill 

U.S.  Army  Vehicle  Propulsion  Directorate/ARL 
Lewis  Research  Center,  Cleveland,  Ohio  44135 


Invited  Presentation 


DOE/NASA  Engine  Research  Program 

M.  Valeo 

NASA  Lewis  Research  Center 
Vehicle  Propulsion  Directorate,  Army  Research  Laboratory 
Cleveland,  OH  44135-3127 
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Invited  Presentation 


Requirements  for  the  Next  Generation  Diesel  Engine 


J.M.  Clarke 

Engine  Research  Department 
Caterpillar  Inc. 

Peoria,  IL  61629 


The  common  goals  of  the  ARO,  universities  and  the  heavy  duty  engine  industry  are  discussed.  It 
is  apparent  that  understanding  emission  mechanisms,  evolving  the  next  generation  of  design 
optimization  methods  and  identifying  potential  revolutions  in  engine  design  should  be  highlights  of 

engine  related  research.  . 

The  chemical  kinetic  and  fluid  dynamics  modeling  of  heterogeneous  combustion  processes 
remains  beyond  the  reach  of  engine  designers  and  developers.  The  low  speed,  uncertain  accuracy 
and  high  cost  of  current  physically  based  simulations  takes  them  outside  the  design  and 
development  cycles.  This  is  the  area  which  needs  the  most  development  of  physical  insight  leading 

to  efficiently  formulated  models.  . 

Current  engine  design  is  an  extremely  complex  decision  process  based  on  information  from 
marketing,  raw  materials  costs,  manufacturing  costs,  and  analyses  related  to  durability  and 
reliability,  performance  and  emissions.  Much  of  the  data  is  subject  to  uncertainty  and  some  of  it  is 
based  on  experience  in  the  form  of  "rules  of  thumb".  In  spite  of  this  informal  process  (some  might 
say  because  of  its  informality)  engines  have  improved  enormously.  Additional  improvement  is 
anticipated  based  on  up-to-date  information  processing  and  rational  optimizing  procedures.  For 
instance,  we  are  a  long  way  from  the  genetic  programming  of  optimum  engine  designs. 
Anticipated  computer  performance  increase  seems  to  bring  procedures  of  this  type  into  the  realm  of 

possibility.  Efforts  to  move  in  this  direction  will  continue  to  attract  industry  interest 

Since  the  diesel  engine  matured  in  its  current  form  there  have  been  many  "breakthroughs 
in  technical  capability.  Looking  at  the  "breakthrough"  list  it  becomes  tempting  to  try  new  engine 
designs  based  on  recognizing  the  role  such  capability  might  play  in  radical  approaches  to  a  new 
generation  of  ecologically  benign  and  economically  attractive  engines.  Here  is  a  chance  for 
universities  to  challenge  their  multi-disciplined(?)  eager(?)  young  engineers.  A  table  of  required 
and  desired  functional  features  in  any  new  design  is  included. 


4-1 


High  Pressure  Preignition  Chemistry  of  Hydrocarbons  and  Hydrocarbon  Mixtures 


N.P.  Cemansky  and  D.L.  Miller 
Department  of  Mechanical  Engineering  and  Mechanics 
Drexel  University 
Philadelphia,  PA  19104 


An  experimental  program  to  study  the  effects  of  pressure  on  preignition  chemistry  of 
hydrocarbon  fuels  is  currently  underway.  Experiments  are  being  conducted  over  a  range  of 
operating  conditions,  including  temperatures  which  encompass  the  low  and  intermediate 
temperature  regimes  (500-1000  K),  and  pressures  ranging  from  less  than  1  atm  to  20  atm.  Three 
experimental  test  facilities  are  being  used  to  perform  this  study:  a  pressurized  flow  reactor  (PFR), 
an  atmospheric  flow  reactor  (APFR),  and  a  static  reactor  (SR).  Results  are  being  used  to  provide 
much  needed  kinetic  and  mechanistic  information  in  these  regimes  over  a  range  of  pressures. 
Additionally  an  effort  to  develop  in  situ  diagnostic  capabilities  in  the  laboratory  has  received 
partial  funding  under  this  contract 

Recent  efforts  have  concentrated  on  the  SR,  PFR,  diagnostic  development  and  modeling. 
Studies  of  n-pentane  in  the  SR  have  identified  important  steps  which  control  product  formation. 
Results  on  1-pentene  indicate  that  long  alkenes  exhibit  negative  temperature  coefficient  behavior 
similar  to  alkanes.  Experimental  studies  of  n-butane,  iso-butane  and  n-pentane  oxidation  at 
elevated  pressures  using  the  PFR  have  confirmed  that  branched  chains  are  less  reactive  than 
straight  chains  and  have  provided  species  distribution  data  for  use  in  modeling.  The  newly 
installed  on-line  FTIR  measured  significant  levels  of  formic  acid,  a  species  difficult  to  measure 
by  other  techniques,  from  the  reaction  of  n-butane,  iso-butane  and  n-pentane.  DFWM 
development  continues  with  measurements  of  OH  and  continued  work  on  H02-  All  of  these  data 
are  used  to  develop  detailed  chemical  mechanistic  information  and  elucidate  the  effects  of 
pressure  on  the  mechanisms.  We  have  recently  adopted  a  two  temperature  sensitivity  method  to 
test  for  the  effect  of  reaction  rate  changes  on  NTC  behavior. 


Research  supported  by  the  Army  Research  Office  Under 
Contract  No.  DAAH04-93-G-0042;  Proposal  No.  30782-EG 
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Data  needed  to  extend  and  validate  models 


Effects  of  Temperature  and  Pressure 
on  Overall  Kinetic  Mechanism _ 
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H  +  02  turnover  marks  the  transition  between  the  intermediate  and  high 
temperature  regimes 


The  Low  and  Intermediate 

Temperature  Regimes  in  Diesel  Engines _ 

Adiabatic  (T,P)-Histories  in  Internal  Combustion  Engines  Showing  the  Importance  of  the 
Low  and  Intermediate  Temperature  Regimes  and  Comparing  the  Experimental 
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-  Development  of  optical  diagnostics  for  on-line  and  in  situ  analysis 

-  Chemical  kinetic  modelling 


Static  Reactor  Activities 

-  Experimental  mechanistic  studies  with  ethene,  propene,  propene  oxide. 


S 

2 

C fl 


C/2 

o 


o 


G 

43 

G 


I 


N 

.G5 

PH 

33 

g 

•g 

S 

■s 

C/2 

<D 


<L> 


G 

<D 

CX 

i 

G 

33 

§ 

C3 

G 


G 

o 

C/2  ^ 


W) 

c 

13 

33 

o 

S 

o 

"X3 

G 


G 


G 

a 

o 

J-h 

G 

rv 

C/5 

<L> 

G 
<D 

G 

#v 

C/2 
0 

•  •  e3 

O  G 

■e  g 

g  -g 

2^ 

33  33 

>*  H 

45  43 

2  C 

u  s 

ti.^ 

G 

Q4h  a 

ZI  W)  g 

<S  'c3  R 

iG  1  2* 

•§  «  s 

3  w  o 

Cfl  lH  O 

c  '3  u 

.2  g  8 

■W  X2 

*a^  ~ 

w)  <u  *3 

•*-i  3  G 
75  «  G 

-a  S  | 


C/2 

2 

S 

x 


G 

0 


C/2 

0 


ci 

0 

& 

£ 

O 

E 

2 

S 

C/2 

0 


•G  G 


•?  £ 

*53  2 

O  *+H 

<  G 


33 

E 

w 

H 

*G 

(D 

•g 

C/5 

•4— >  0 

c/2  C 

§ 

i 

0 

0 

43 

a 

§ 

rl 

r-4 

0  s 

>  4-J 

U  J 


C/2 

O 

s 


0  £  £8 


C/2 

0 


O 

< 

i-t 

o 

g 

G 

0 

P4 

£ 

E 

<L> 

N 

'§ 

C/2 

C/2 

2 

cm 


5-6 


-  Reactivity  mapping  of  C2-C4  alkanes  and  alkenes 

-  Detailed  species  measurements  of  propane  and  n-butane  oxidation 

-  Examined  the  effect  of  pressure  on  the  turnover  between  the  low  and 
intermediate  temperature  regimes 
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Static  Reactor  Recent  Work 
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Further  analysis  will  provide  insight 
on  NTC  region  transition  in  olefins 
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Iso-butane  oxidation  studies 

0  McCormick,  M.S.  thesis,  June  1994 

n-Pentane  oxidation  studies 

0  Wood,  M.S.  thesis,  August  1994 
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-  Optimizing  beam  tilting  technique  to  increase  sensitivity 

-  Working  on  CN  and  H02  measurements 

-  Developing  a  technique  for  H02  calibration  using  electrical 
discharge  under  low  pressure  in  a  discharge  flow  reactor 


Chemical  Kinetic  Modeling 

•  Current  detailed  kinetic  models  have  difficulty  reproducing 
experimental  data  that  span  multiple  reaction  regimes  (e.g.,  the  NTC  region) 

•  In  cooperation  with  LLNL  we  are  working  to  improve  these  models  by: 
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Status  and  Issues 
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Measurement  of  key  radical  intermediates  may  simplify  the  puzzle 
-  DFWM  has  reasonable  potential  for  measuring  H02  radicals 


Current  and  Planned  Activities 
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-  Continue  mechanistic  development  activity 

-  Continue  collaborative  efforts  with  LLNL 


Ignition  in  Convective-Diffusive  Systems 


C.  K.  Law  and  T.  G.  Kreutz 
Department  of  Mechanical  and  Aerospace  Engineering 
Princeton  University 
Princeton,  NJ  08544 


Autoignition  in  non-prenuxed  hydrocarbon/air  combustion  has  been  experimentally,  numerically, 
and  analytically  investigated  in  counterflowing  jets  of  fuel  against  hot  air.  Details  of  die  ignition 
process  are  studied  by  examining  the  steady-state  situation  as  the  air  temperature  is  raised  to  a  point 
just  prior  to  ignition.  Experimental  measurements  of  the  ignition  temperature  have  been  made  as  a 
function  of  system  pressure,  aerodynamic  strain  rate,  and  fuel  concentration  for  fuels  such  as 
hydrogen,  methane,  propane,  n-butane,  and  iso-butane.  These  studies  are  augmented  by  numerical 
modelling  using  detailed  chemistry  that  elucidates  the  nature  of  underlying  ignition  ‘kernel’,  a 

spatially-localized  zone  of  high  chemical  reactivity. 

In  the  H2-air  system,  both  experiments  and  calculations  show  that  the  ignition  temperature, 
when  plotted  as  a  function  of  pressure,  exhibits  a  characteristic  ‘Z’-shape  similar  to  that  seen  in  the 
explosion  limits  of  homogeneous  H2-air  mixtures.  The  three  sections  of  this  Z-curve  are  denoted 
the  first,  second,  and  third  ignition  limits  in  analogy  with  the  three  explosion  limits.  Numerical 
analysis  indicates  that,  for  each  limit,  the  dominant  chemistry  is  identical  in  both  the  explosion  and 
ignition  situations.  Aerodynamic  straining  is  found  in  all  cases  to  increase  the  ignition  temperature 
due  to  the  decreased  size  of  the  ignition  kernel  and  the  decreased  residence  time  of  active  radicals 
within  it  Ignition  in  the  first  and  third  limits  are  found  to  be  sensitive  to  aerodynamic  straining 
because  loss  of  radical  species  via  mass  transport  plays  an  important  role  in  these  two  ignition 
regimes.  In  contrast,  the  second  ignition  limit  is  governed  primarily  by  fast  chemical  reactions  and 
is  thus  relatively  insensitive  to  variations  in  the  strain  rate. 

The  chemistry  governing  the  H2-air  ignition  process  has  been  distilled  by  means  of  careful 
analysis  and  verified  by  means  of  ‘skeletal’  chemical  reaction  mechanisms  consisting  of  a  minimal 
number  of  elementary  reaction  steps.  The  skeletal  mechanisms  have  been  further  compressed  to 
‘reduced’  mechanisms,  consisting  of  non-elementary  reactions,  via  techniques  such  as  the  steady- 
state  approximation.  These  methods  and  reaction  schemes  are  validated  by  means  of  comparison 

with  the  results  of  computations  using  the  full  kinetic  mechanism. 

In  the  H2-air  system,  the  role  of  heat  release  in  providing  ‘thermal  feedback’  at  the  ignition 
turning  point  has  been  examined  in  detail  for  all  three  ignition  limits.  Contrary  to  classical  notions 
based  upon  one-step  overall  chemistry,  thermal  feedback  is  shown  to  play  essentially  no,  or 
minimal,  role  in  the  steady-state  solution  at  the  ignition  turning  point  -  either  in  its  character  or 
parametric  dependence.  In  most  cases,  turning  point  and  S -curve  behavior  are  found  to  exist  in  the 
complete  absence  of  heat  release,  driven  solely  by  ‘kinetic’  feedback  provided  by  algebraic 
nonlinearities  in  the  coupled  chemical  kinetics.  As  a  result,  the  location  of  the  ignition  turning 
points  are  found  to  be  essentially  governed  by  the  kinetics  of  gain  versus  loss  of  key  radicals  in  the 
ignition  kernel  which  depend  parametrically  upon  global  variables  such  as  air  temperature,  strain 

rate,  pressure,  and  fuel  concentration.  . . 

Experimental  and  calculated  ignition  temperatures  for  methane,  as  well  as  measured  ignition 
temperatures  for  n-butane  and  iso-butane,  are  found  to  be  quite  sensitive  to  the  strain  rate. 
Sensitivity  to  aerodynamic  straining  is  found  to  occur  when  at  least  one  critical  ignition  reaction 
proceeds  at  a  rate  which  is  relatively  slow  and  thus  comparable  to  the  rates  of  mass  transport  out  of 
the  localized  region  of  the  ignition  kernel. 
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Background 


•  Ignition  in  Diesel  engines  is  expected  to  involve 
inhomogeneous  mixtures: 

-  Premixed  reactants 

-  Nonpremixed  reactants 

-  Temperature  nonuniformity 

•  Presence  of  temperature  and  concentration  gradients 
implies  the  importance  of  diffusive  transport. 


•  Ignition  states  determined  from  homogeneous 
experiments  may  not  be  directly  applicable. 


Objectives 


To  study  the  dynamics  and  chemical  kinetics  of 
hydrocarbon/air  ignition  in  the  presence  of  diffusive 
mixing. 

Experimentally  and  computationally  determine: 

1)  Bulk  ignition  parameters  such  as  ignition 
temperature  and  strain  rate; 

2)  Identify  the  critical  kinetic  steps  and  species 
governing  diffusive  ignition. 

Analytically  derive: 

1)  Reduced  mechanisms  relevant  to  diffusive  ignition; 

2)  Ignition  Damkohler  number  and  criteria. 

Conceptually  understand/unify: 

1)  Homogeneous  versus  diffusive  ignition; 

2)  Thermal  versus  kinetic  ignition. 


Approach 


•  Counterflow  configuration:  well-defined  strain  rate  and 
hence  characteristic  flow  time 

•  Steady-state  ignitabiiity  studies 

•  Transient  ignition 

•  Nonpremixed  ignition:  cold  fuel  against  hot  air 

•  Premixed  ignition:  cold  fuel/air  against  hot  inert 

•  Investigate/manipulate  ignitabiiity  for  different: 

-  Fuels  (e.g.  H2,  CH4,  C3H8,  n-C^io,  FC4H-10) 

-  Fuel  concentration 

-  Temperature  of  the  hot  jet 

-  System  pressure 

-  Flow  strain  rate 


Diffusive  Ignition  of 
Hydrogen/Oxygen  Mixtures 


Homogeneous  ignition  of  hydrogen/oxygen  mixture  is 
characterized  by  three  explosion  limits. 

How  is  this  behavior  modified  by  diffusion  and 
aerodynamic  straining  in  a  nonuniform,  flowing  mixture? 


H 2/Air  Reaction  Mechanism 

(Yetter,  Dryer,  Rabitz,  1991) 


Reaction 

A 

n 

Ea 

1. 

H  +  02  <=>  0  +  OH 

1.92x1014 

0.0 

16.44 

2. 

0  +  H2  «  H  +  OH 

5.08x1 04 

2.67 

6.29 

3. 

OH  +  H2  «  H  +  H20 

2.16x10® 

1.51 

3.43 

4. 

OH  +  OH  <=>  0  +  H20 

1.23x104 

2.62 

-1.88 

5.* 

H2  +  M  <=>  H  +  H  +  M 

4.57x1 019 

-1.4 

104.4 

6.* 

O  +  O  +  M  <=>  O2  +  M 

6.1 7x1 015 

-0.5 

0.0 

7.* 

O  +  H  +  M  oOH  +  M 

4.72x1 018 

-1.0 

0.0 

8.* 

H  +  OH  +  M  <=>  H2O  +  M 

2.25x1 022 

-2.0 

0.0 

9.* 

H  +  O2  +  M  <=>  HO2  +  M 

6.1 7x1 019 

-1.42 

0.0 

10. 

HO2  +  H  <=>  H2  +  O2 

6.63x1 013 

0.0 

2.13 

11. 

H02  +  H  <=>  OH  +  OH 

1.69x1014 

0.0 

0.87 

12. 

HO2  +  0  0  OH  +  O2 

1.81x1013 

0.0 

-0.4 

13. 

HO2  +  OH  <=>  H2O  +  O2 

1 .45x1 016 

-1.0 

0.0 

14. 

HO2  +  HO2  <==>  H2O2  +  O2 

3.02x1 012 

0.0 

1.39 

15.* 

H2O2  +  M  <=>  OH  +  OH  +  M 

1 .20x1 017 

0.0 

45.5 

16. 

H2O2  +  H  <=>  H2O  +  OH 

1 .00x1 013 

0.0 

3.59 

17. 

H2O2  +  H  <=>  H2  +  HO2 

4.82x1 013 

0.0 

7.95 

18. 

H202  +  O  0  OH  +  H02 

9.55x10® 

2.0 

3.97 

19. 

H2O2  +  OH  <=>  H2O  +  HO2 

7.00x1 012 

0.0 

1.43 

Third  body  enhancement  factors:  H20: 12;  H2: 25 
Units:  cm,  mole,  sec,  kcal 
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Peak  Mole  Fraction 
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6-9 


Pressure  (atm) 
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Thermocouple 

indicator 


Ignition  Temperature 
(Radiation  Corrected) 
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Reduced  Mechanisms  for 
Hydrogen/Oxygen  Ignition 

•  First  Limit:  O  and  OH  in  steady  state 

3H2  +  O2  ” ^  2H2O  +  2H 
H  +  O2  +  M  — >  HO2  +  M 
2HO2  — ^  H2  +  2O2 

•  Second  Limit:  H,  O,  and  OH  in  steady  state;  2coa  ~  cob 

H  +  02  ->  O  +  OH  (a) 

H  +  O2  +  M  — >  HO2  +  M  (b) 


•  Third  Limit:  H,  O,  OH,  and  H02  in  steady  state 

3H2  +  2O2  — ^  2H2O  +  H2O2 
2H2  +  O2  2H2O 
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Methane/Air  lanition:  i-Butane/Air  Ignition 
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CFL  Mole  Fraction  Mole  Fraction  iso-Butane 


Summary  of  Accomplishments 


•  Nonpremixed  ignition  of  hydrogen/air  extensively 
studied; 

•  Unified  chain-transport  interpretation  of  ignition 
phenomena; 

•  Role  of  thermal  versus  kinetic  ignition  identified; 

•  New  numerical  algorithm  developed  for  the  study  of 
ignition/extinction  phenomena. 


Current  and  Future  Work 

•  Premixed  systems 

•  Transient  ignition 

•  Methane  and  methane/hydrogen  ignition 

•  C4-hydrocarbon  ignition 
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Resonant  Holographic  Interferometry 


P.A.  DeBarber,  N.J.  Brock,  M.S.  Brown,  C.F.  Hess,  and  J.D.  Trolinger 

MetroLaser 

18006  Skypark  Circle,  Suite  108 
Irvine,  CA  92714 


This  presentation  discusses  a  novel,  nonintrusive,  spectroscopic  tool  to  explore 
combustion  and  flames.  This  tool  is  called  Resonant  Holographic  Interferometric 
Spectroscopy  (RHIS).  RHIS  is  a  completely  nonintrusive  optical  diagnostic.  As  the 
RHIS  acronym  implies,  RHIS  is  a  hybrid  technique.  It  brings  together  the  three- 
dimensional  imaging  capability  of  holography ,  the  phase  sensitive  detection  of 
interferometry ,  and  the  species  specificity  of  spectroscopy  to  form  a  powerful,  versatile 
diagnostic  package.  RHIS  is  a  holographic  interferometric  technique  which  may  be 
combined  with  tomography  to  provide  three  dimensional  characterization  and 
quantitative  measurement  to  species  number  density,  temperature,  velocity,  pressure, 
molecular  mixing,  and  other  thermophysical  quantities. 

Conceptually,  a  single  RHIS  measurement  records  the  concentration  of  a 
resonant  absorbing  chemical  species  in  the  form  of  a  spatially  resolved  interferogram. 
This  is  achieved  by  holographically  recording  the  phase  shift  due  to  the  resonant 
species  simultaneously  with  the  phase  shift  due  to  all  other  nonresonant  species.  Upon 
reconstruction,  the  phase  shift  produced  nonresonantly  is  subtracted  out  leaving  a  neat 
interferogram  which  contains  information  only  about  the  resonant  species.  In  practice, 
this  is  achieved  in  the  following  way.  The  output  beams  of  two  lasers  operating  at 
different  frequencies  are  spatially  and  temporally  overlapped  to  form  a  single  two-color 
beam.  One  of  the  lasers  is  tuned  to  a  molecular  resonance  of  the  species  of  interest. 
The  other  laser  is  tuned  off-resonance.  The  two-color  beam  is  then  split  into  the 
standard  object  and  reference  legs  of  a  conventional  holographic  interferometer.  Light 
from  each  leg  is  then  recombined  on  a  film  plate  to  form  two  holograms  simultaneously 
at  each  frequency.  Both  of  the  two  electric  fields  at  the  two  different  frequencies  in  the 
object  arm  are  phase  shifted  by  the  same  amount  due  to  density  fluctuations.  However, 
the  resonantly-tuned  electric  field  receives  an  additional  phase  shift.  Upon 
reconstruction,  the  two  holographically  recorded  images  interfere  with  each  other 
resulting  in  a  cancellation  of  all  phase  information  that  was  common  to  the  two  electric 
fields.  Therefore  resulting  in  an  interferogram  which  preserves  only  the  phase  shift 
fringes  attributable  to  the  resonant  species.  This,  in  effect,  is  a  background-free 
recording  of  the  spatially  resolved  density  of  the  resonant  species.  Since  phase  shifts 
are  largest  near  resonance,  the  RHIS  technique  is  far  more  sensitive  to  low  species 
concentrations  than  conventional  nonresonant  interferometry.  The  detected  phase  shifts 
are  linearly  dependent  on  concentration  which  makes  the  RHIS  technique  also  linearly 
dependent  on  concentration. 


7-1 


on  resonance 


Figure  1.  Plot  of  Real  and  Imaginary  components  of  the  index  of  refraction  near  an 
absorption  line. 


Figure  2.  Two  wavelength  RHIS  procedure.  Panel  (a)  illustrates  the  recording  procedure. 

corresponds  to  the  laser  tuned  "on  resonance";  X2  corresponds  to  the  laser  tuned  "off 
resonance".  Panel  (b)  outlines  the  reconstruction  process.  Note  that  the  reconstruction  can 
be  accomplished  with  almost  any  available  laser  >4. 
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Figure  3.  Schematic  diagram  of  the  two-laser  RHIS  system.  In  this  case,  the  object 
probed  is  the  flame  produced  in  a  laboratory  burner. 


B.  Conventional 
Interferogram 


Figure  4.  (a)  Experimental  configuration,  (b)  Conventional  holographic  interferogram. 
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A.  Two-wavelength,  off- 
resonance  Holographic 
Interferogram 


B.  RHIS  Holographic 
Interferogram  of  Na 


Figure  5.  (a)  two-wavelength  off-resonance  interferogram  illustrating  total  subtraction 
inherent  in  holographic  reconstruction,  (b)  RHIS  interferogram  of  Na  seeded  into  a 
methane/air  diffusion  burner  flame. 


Figure  6.  Predictions  for  the  dynamic  range  of  an  individual  RHIS  hologram  as  a  function 
of  the  minimum  detectable  fringe  shift. 
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CONTROL  OF  SOOT  FORMATION 

R.A.  Dobbins 
Division  of  Engineering 
Brown  University 
Providence,  RI  02912 

Micrographs  (TEM)  obtained  by  thermophoretic  sampling  of  soot  particles  from  diffusion 
flames  have  revealed  two  contrasting  morphologies.  On  the  fuel  side  of  the  lower  diffusion 
flame  front,  many  small  singlet  precursor  particles  are  found  that  are  nearly  transparent  to  the 
electron  beam.  In  the  higher  temperature  regions  the  particles  are  found  to  consist  of  opaque, 
carbonaceous  aggregates  made  up  of  monodisperse  spherules.  Our  studies  and  reference  to  the 
literature  lead  to  the  conclusion  that  the  precursor  stage  observed  in  diffusion  flames  is  a 
common,  if  not  universal,  stage  in  the  formation  of  the  carbonaceous  materials  in  shock  tubes,  in 
premixed  and  diffusion  flames,  and  in  droplets  in  micro  gravity. 

The  nature  of  the  precursor  particle  recently  has  been  explored  by  laser  microprobe  mass 
spectrometry  (LMMS).  This  research  has  been  a  collaboration  with  R.A.  Fletcher  of  the  NIST 
using  the  LAMMA-500  instrument  available  there.  The  advantages  and  drawbacks  of  the 
LMMS  method  and  the  results  of  testing  samples  of  pure  and  mixed  polycyclic  aromatic 
hydrocarbons  are  reviewed.  LMMS  mass  spectra  of  precursor  material  from  low  on  the  axis  of 
the  diffusion  flame  reveals  many  ion  masses  corresponding  to  PAHs  in  the  200  to  300  amu 
range.  The  most  prominent  of  these  are  the  mass  numbers  are  252,  e.g.,  benzo(a)pyrene;  276, 
e.g.,  benzo(ghi)perylene;  and  300,  coronene.  Of  striking  significance  is  that  these  PAH  masses 
were  predicted  by  Stein  and  Fahr  (1985)  to  be  among  the  most  thermodynamically  stable  species 
("stabilimers")  in  typical  flame  conditions  of  pressure  and  temperature.  Mass  spectra  of  the 
carbonaceous  material  from  the  upper  flames  show  many  carbon  clusters  with  typically  with  0, 
1,  or  2  hydrogen  atoms.  Small  quantities  of  ion  masses  in  the  418  to  444  amu  range  are  also 
found  in  the  carbonaceous  materials  along  with  traces  of  fullerene-like  ions  -  600, 720,  840  and 
1105  amu.  The  values  of  hydrogen  mole  fraction,  or  C/H  ratio,  of  both  precursor  and 
carbonaceous  material  are  found  by  the  LMMS  to  conform  to  prior  expectations. 

The  rate  of  carbonization  of  precursor  materials  is  an  important  parameter  in  the 
understanding  and  intervention  in  the  process  of  soot  formation.  In  nitrogen-diluted  ethene 
diffusion  flames  the  onset  of  carbonization  is  delayed  and  occurs  over  a  range  of  temperatures 
according  to  the  amount  of  dilution.  Thermocouple  measurements  are  conducted  using  the 
rapid  insertion  technique  to  eliminate  the  influence  of  soot  deposition  on  sensor  emissivity.  The 
observed  locations  of  carbonization  and  the  observed  temperature  profiles  lead  to  carbonization 
rate  measurements  when  first  order  kinetics  are  assumed.  Two  methods  of  data  reduction  are 
employed  and  give  similar  results.  An  activation  energy  of  21  kcal/mole  is  found  for  the  rate  of 
carbonization  of  precursor  particles.  The  implications  of  this  carbonization  rate  on  soot 
formation  in  combustion  processes  are  discussed. 


ARO  Grant  DAAL03-92-G-0023 
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Evolution  of  Soot  in  Flames 


LASER  MICROPROBE  MASS  SPECTROMETRY 


-  265  NM  LASER  OUTPUT  FOCUSSED  ON  SAMPLE  - 1  TO  100  MICROJOULES, 

10ns  PULSE. 

-  BEAM  DEFOCUSSED  TO  50-  100  MICRON  DIAMETER,  10+8  W/m2. 

-VACUUM  OF  10-4  PA. 

-  MATERIAL  IS  SIMULTANEOUSLY  ABLATED  AND  IONIZED. 

-  IONIZED  PRODUCTS  MASS  ANALYZED,  0  TO  5700  AMU,  M/AM  -500  @  M/Z  = 

200,  M=1  TO  5700  AMU. 

-  REPRODUCIBILITY  LOW  BECAUSE  MICROSCOPIC  VARIABILITY  OF  SAMPLE 

COMPOSITION  AND  THICKNESS  &  VARIABILITY  OF  LASER  PULSE. 
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LMMS  Analysis  -  Soot  from  C2H4  diffusion  flame. 

Z  =  50  mm  T  =  1570  K 

x  ®  4  to  24 
ManyCxHy  y  =  0,1,2 

Xh  >  0.14  C/H  *  6.1 

Small  amounts  of  PAH-like  species  418-444 
Trace  levels  of  masses  corresponding  to  C50,  C60,  C70,  C92 

Z  =  20  mm  T=  1430  K 

PAH-like  ions  (count  fraction) 

252  (0.16)  276  (0.20)  300  (.090) 

Xh^0.35  C/H  «  1.9 

Many  CxHy  Clusters 

Ten  PAHs-like  ions  reported  by  others 
Five  PAHs-like  ions  not  reported 
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Species 

Formula 

M 

(g/mole) 

Cyclopen  ta(cd)pyrene 

Cl8H]0 

226 

Benzo<a)pyrene 

C20H12 

252 

Benzo(gh/)perylene 

C22H12 

276 

Indeno(l,2,3-af)pyrene 

C22HI2 

276 

Anthanthrene 

C22H12 

276 

Coronene 

C24H12 

300 

Naphtho(8, 1 ,2-<2Ac)coronene 

C30HI4 

374 

Ovalene 

C32HU 

398 

no 


no 


goo 


jOOD 

DOJ 


OQ 

□or 


jOOOr 

Ow0 

lOOOr 


Table  3.6-2:  PAH  species  identified  in  flame  samples  using  liquid 
chromatography.  Pyrene  and  fluoranthene,  which  were  also  detected  by 
gas  chromatography,  are  not  listed  in  this  table.  Many  peaks  which 
appeared  in  the  chromatographs  were  not  identified. 

(From  Feitelberg,  Ph.D.  Thesis,  MIT,  1990) 

PAHs  Found  in  C2H4  Premixed  Flame 
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Ethene  Diluted  with  Nitrogen  (Q  of  C2H4  =  2.1  cm3/s) 


2000 


Figure  1.  Diffusion  Flame  Axial  Temperature  Profiles 
for  Pure  and  Diluted  Ethene 


CARBONIZATION  RATE 


METHOD  I 

Kca  =  [t(Z2)  -•(ZpJ1 

Zi,  Z2,  observed,  5  mm  resol. 
T(z)  for  Six  Dilution  Ratios 

Plot  ln(Kca)  vs.  T'1,  find  A  &  E 


METHOD  H 

-  A,  '  "  ' 


Given  ti,  T(t),  assume  E,  reiterate  to  find  A  to  give 


=  e-l  when  t 
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Soot  Modeling  and  Experiments  at  the  ERC 

D.E.  Foster 

Engine  Research  Center 
University  of  Wisconsin  -  Madison 

In  this  presentation  the  Engine  Research  Center's  experimental  and  modeling  work  relative 
to  soot  emissions  from  diesel  engines  will  be  reviewed.  First  an  overall  picture  of  the 
work  wiU  be  presented  which  wiU  be  foUowed  by  a  focused  discussion  on  three  specific 
projects.  The  three  projects  being:  in-situ  measurements  of  soot  particle  size,  volume 
fraction  and  number  density  via  a  combination  two  wavelength  scattering/extinction 
technique,  the  development  of  a  phenomenological  soot  model  for  the  purpose  of 
inclusion  into  three  dimensional  computational  diesel  simulations  and  the  effort  to  predict 
soot  emission  with  simplistic  soot  models  taken  from  the  literature  and  applied  to  the 

KTVA II  simulation.  . 

The  scattering  extinction  technique  was  built  on  our  experience  of  installi  g 

windows  into  an  operating  diesel  engine  by  introducing  a  laser,  via  fiber  optic 
transmission,  into  the  cylinder  next  to  the  window.  The  technique  was  compared  to  a 
scattering-two  color  radiant  emission  technique  and  the  results  for  particle  size,  number 
density  and  volume  fraction  were  found  to  be  consistent.  The  measurement  system  was 
then  applied  to  the  engine  operating  under  different  conditions. 

The  phenomenological  soot  model  describes  the  processes  of  inception,  nucleation, 
surface  growth,  oxidation  and  coagulation  with  global  kinetic  reactions.  The  rate 
constants  for  each  global  reaction  is  taken  to  be  consistent  with  those  for  similar  types  of 
reactions  from  published  data.  The  model  has  been  applied  for  conditions  of  a  laminar 
flame  published  in  the  literature  and  the  results  were  encouraging.  Next  the  model  has 
been  exercised  for  a  computational  cell  from  a  KIVA  diesel  engine  simulation.  e 
profiles  of  fuel,  oxygen  and  temperature  were  scaled  in  a  parametric  study  to  examine  the 
response  of  the  model.  No  parameters  were  adjusted  during  these  variations.  The  mo  e 

predicts  the  expected  soot  emission  behavior. 

Finally  the  results  of  our  efforts  to  incorporate  the  Hiroyasu  soot  model  into  the 
KTVA  II  diesel  engine  will  be  discussed.  Hiroyasu's  model  has  been  modified  to  include 

the  Nagle  Strickland-Constable  oxidation.  The  results  have  been  tuned  for  a  smgle 

operating  condition  and  the  model  used  to  assess  whether  it  predicts  the  correct  trends 
with  changes  in  engine  operating  conditions.  The  results  reinforce  the  strong  coup  g 
between  the  combustion  model  and  the  emission  models.  It  has  been  determined  that  the 
NOx  model  demonstrates  the  highest  level  of  sensitivity  to  small  variations  in  the 

combustion  model. 
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"Soot  Modeling  and  Experiments  at  the  ERC" 


Dave  Foster 
Engine  Reserach  Center 
University  of  Wisconsin  -  Madison 

Eleventh  ARO  Engine  Workshop 
March  15, 1994 

ERC  Faculty:  Gary  Borman,  Mike  Corradini,  Pat  Farrell,  Dave  Foster,  Jay 

Martin,  Rolf  Reitz  and  Chris  Rutland 

Overview  of  Current  Work: 


Experiments: 

.  Two  dimensional  two  color  optical  pyrometry 
In-situ  measurement  of  particle  size,  number 
density,  volume  fraction  and  temperature  via 
scattering  and  extinction 
.  Interaction  between  multiple,  split  injection 
scenarios  and  the  soot  NOx  trade-off 


Modeling  Efforts: 

.  Development  of  a  phenomenological  soot  model 
.  Incorporation  of  literature  models  into  KI V A 

.  Implimentation  of  the  more  detailed  models  into 
KIVA 


In-situ  Measurement  of  Particle  Size,  Number  Density, 
Volume  Fraction  and  Temperature  via  Scattering  and 

Extinction 


.  Point  measurement  in  the  cylinder  of  an  operating  DI 
diesel  engine 

.  Simultaneous  measurement  of  light 

scattering/extinction  and  radiant  emission  from  engine 

soot  cloud 

.  Different  engine  operating  conditins  were  investigated 


Phenomenological  soot  model 

.  Global  kinetic  expressions  for  the  physical  processes 
of: 

inception 
surface  growth 
oxidation 
coagulation 

•  Rate  constants  for  the  global  kinetic  expressions  are 
consistent  with  literature 


Compared  to  laminar  flame  data 
Implimented  in  a  KIVA  computational  cell 


Incorporation  of  literature  models  into  KIVA 

Hiroyasu  model  with  NSC  oxidation  implimented  into 
KIVA  II  diesel  engine  simulation 

Preparations  for  more  detailed  models  to  be 
implimented  into  the  simulation 


Schematic  diagram  of  optical  set-up  and  engine 


0.508mm  O.D. 
Surgical  Tubing 

Copper  Gasket 


3mm  Spherical 
Glass  Lens 


Screw 


SENDING  WINDOW 
UNIT 


1.27  X  1.27mm  Slot 


1.02mm  O.D. 
Surgical  Tubing 


Copper 

Gaskets 


Sender  Window  Hole 


Copper  Gasket 


Retaining 

Nut 


Random  Fiber 
Optic  Bundle 


Aperture 

Lenses 

INVAR 

Sapphire 

Window 

Crevice 

Volume 


Cross  section  of  the  optical  probe. 


Cross  section  of  the  engine  head,  piston  bowl,  and  measurement  volume 
location. 
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d-q EH3 □  Radiation  850  nm 
a-a-a-a-a  Radiation  514  nm 
00000  Scattering  514.5  nm 
ii:i)  Scattering  488  nm 


Separated  radiation  and  scattering  signals  after 
ensemble  averaging  at  the  three  measured  wavelengths. 
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Temperature  (K)  Diameter  (nm) 


Soot  Volume  Fraction,  fv 


Soot  Volume  Fraction,  fv 


concentration  [mol/cm3] 


inert 

products 

f 

(3)  1  oxidation 
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REGIMES  OF  JET  BREAKUP 
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ABSTRACT 

The  regimes  of  jet  breakup  including  the  first-wind  induced,  second  wind-induced,  and 
atomization  regimes  are  elucidated  by  use  of  a  single  equation.  The  equation  is  obtained  by 
use  of  linear  stability  analysis.  The  equation  relates  the  dimensionless  intact  length  L  to  the 
dimensionless  velocity  U  for  relevant  flow  parameters.  The  results  are  plotted  in  the  L  vs  U 
space  to  elucidate  the  regimes  and  mechanism  of  the  jet  breakup.  The  theory  compares  well 
with  some  experiments. 

INTRODUCTION 

The  atomization  process  is  widely  used  in  the  formation  of  sprays  for  various  industrial 
applications,  including  fuel  spray.  The  process  can  be  achieved  by  breaking  up  a  liquid  jet 
or  a  liquid  sheet  emanating  from  a  nozzle  of  various  geometries.  The  main  cause  of 
atomization  is  known  to  be  due  to  the  pressure  fluctuation  [1,2].  The  jet  may  also  be  broken 
up  by  the  Rayleigh  mode  of  capillary  pinching  which  produces  droplets  of  diameter  much 
larger  than  that  of  atomized  droplets  [3].  For  a  given  jet  breaking  up  in  a  given  gas,  the 
Rayleigh  mode  occurs  at  low  jet  speed  and  the  atomization  process  takes  place  at  relatively 
large  velocity.  In  the  Rayleigh  mode  of  breakup,  the  actual  breakup  does  not  take  place  until 
some  finite  distance  down  stream  of  the  nozzle.  This  distance  is  called  intact  length.  The 
intact  length  increase  with  the  jet  speed  up  to  a  critical  speed,  and  then  decreases  as  the  speed 
is  increased  further.  This  transition  takes  place  in  so  called  first  wind-induced  breakup  regime 
[4]  which  follows  the  Rayleigh  breakup  regime.  After  the  transition,  the  jet  continues  to 
shorten  and  numerous  small  ligaments  emanate  from  the  liquid-gas  interface  in  the  so  called 
second  wind-induced  breakup  regime.  Finally  the  intact  length  vanishes,  and  the  spray  of  small 
droplets  is  formed.  Attempts  have  been  made  to  explain  the  change  of  the  intact  length  with 
the  jet  speed  in  different  regimes,  based  on  the  first  principle.  Weber  [5]  attempted  to  explain 
the  first  wind-induced  regime  on  the  basis  of  a  simplified  characteristic  equation  obtained  from 
linear  stability  analysis.  His  theory  predicted  a  considerably  smaller  critical  velocity  for  the 
maximum  intact  length.  Stirling  and  Sleicher  [6]  improved  Weber’s  theory  by  use  of  an 
empirical  correlation.  However,  a  theory  which  explains  all  regimes  of  the  jet  breakup  from 
a  single  equation  has  not  yet  appeared.  Here  we  offer  such  a  theory. 
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THEORY 

The  starting  point  of  our  explanation  is  the  following  equation 


,  ..  o  2k2 

(o)  -  ik)2  +  — 

x  (q  -  ik)  + 


( 


l/(k)  2  k  A.  Wlfa) 

I0(k)  "  X2  +  k2  W) 


G)2  Q 


k2) 


X2  -k2  KpCk)!^) 
X2  +  k2  K1(k)I0(k) 
X2  -  k2  1,00 


X2  ♦  k2  Uk) 


»  0, 


(1) 


R  =  Reynolds  number  =  U^o/v, 

We  =  Weber  number  =  p!  U02  Tq/S 
Q  =  Density  ratio  =  p^Pi 
X  =  k2  -  Re  (G3  -  ik) 

Eq.  (1)  is  the  characteristic  equation  obtained  [4]  from  the  stability  analysis  of  a  liquid  jet 
subject  to  disturbances  whose  Fourier  components  have  a  factor  exp  [i  k  y  +  tot],  where  k  = 
k,  +  ilq  is  the  complex  wave  number,  co  =  cor  +  io)i  is  the  complex  wave  frequency,  y  is  axial 
distance  measured  from  the  nozzle  tip.  In  Eq.  (1),  K  and  I  are  the  modified  bessel  functions 
of  orders  indicated  by  their  subscripts,  p  is  the  density,  v  is  the  kinematic  viscosity,  S  is  the 
surface  tension,  U0  is  the  average  jet  velocity,  r0  is  the  jet  radius  and  the  subscripts  1  and  2 
associated  with  the  physical  properties  denote  respectively  the  liquid  and  gas  phases.  Weber 
used  a  drastically  simplified  equation  of  Eq.  (1).  His  results  did  not  enjoy  agreement  with 
experiments.  For  a  given  jet  at  different  U0,  Sterling  and  Sleicher  [6]  multiply  the  Q-term  in 
Eq.  (1)  by  0.175  and  solve  the  equation  numerically  for  [G3r(k)]m„  assuming  lq  =  0.  The  intact 
length  L  is  then  obtained  from  the  relation  L  =  constant/[G3r(k)]m!K.  The  factor  0.175  was  used 
so  that  their  theoretical  prediction  agrees  best  with  the  experiments.  They  argued  that  this 
multiplication  factor  is  necessary  because  the  air  viscosity  which  may  have  a  damping  effect, 
is  neglected  in  the  theory.  They  also  quoted  the  related  work  of  Benjamin  and  stated  that  the 
air  viscosity  effect  is  much  smaller  than  the  inertia  force  effect  of  gas  at  the  interface.  Thus 
there  is  an  apparent  inconsistency.  If  the  air  viscosity  effect  is  mall,  then  the  reduction  of  the 
gas  inertia  force  represented  by  the  Q-term  in  Eq.  (1)  should  not  have  been  reduced  by  a  large 
factor  of  (1  -  0.175). 

Here  we  will  not  multiply  the  Q-term  by  any  empirical  factor.  However  we  shall  use  a 
rational  approximation  to  simplify  the  mathematics,  and  bring  out  the  physical  mechanisms 
of  the  jet  breakup.  Since  both  phenomena  of  First  and  second  wind-induced  regimes  are  due 
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to  the  interfacial  pressure  fluctuation  [1,2],  and  the  viscosities  play  minor  roles  [3],  we  shall 
neglect  the  Re-term  in  Eq.  (1).  This  is  justifiable  for  the  large  Reynolds  number  jet  which  is 
encountered  in  most  of  the  applications.  With  this  approximation,  Eq.  (1)  is  reduced  to 

©2  +  (co  +  ik)2  QA  =  We-‘k(i  -  k2)  B  =  0,  (2) 

where 

a  ^(kjl^k) 

"  Kj(k)I0(k)  ’ 

R  _  m 

w 

We  define  the  intact  length  L  to  be  the  distance  from  the  nozzle  tip  to  a  downstream  position 
where  the  disturbance  amplitude  has  grown  from  the  initial  value  to  an  amplitude  a  for 
breakup.  Then  a  =  ao  exp(kjL)  and  thus 

L  =  kf1  €n  (a/ao) .  (3) 

kj  in  Eq.  (3)  must  be  found  from  Eq.  (2)  as  the  maximum  amplification  rate  for  a  given  Q  and 
We.  The  numerical  task  required  is  very  considerable  since  the  complex  k  appears  in  the 
arguments  of  the  modified  Bessel  functions.  Fortunately,  the  spatial  and  temporal  growth  rates 
are  related  [7]  by 


k,  =  CD,  +  0(We'2)  , 

(4) 

k,  =  CD;  +  0(We'2)  . 

(5) 

Moreover,  both  the  first  and  second  wind-induced  breakup  occur  at  We'1  <  Q,  [1].  For  a  water 
jet  in  air,  Q  =  0.13  at  100  atm.  Thus  we  may  write  Eq.  (3)  as 

L  =  co,'1  fn  (a/ao)  , 

with  an  error  of  less  than  2%  if  the  jet  is  issued  into  air  of  pressure  less  than  100  atm.  cor 
must  be  obtained  from  Eq.  (2)  as  the  maximum  temporal  growth  rate  co^  for  given  Q  and  We. 
The  real  and  imaginary  parts  of  Eq.  (2)  are 

(cor2  -  co,2)  (1  +  QA)  -  kQA(2cO;  +  k)  -  We’1  k(l  -  k2)  B  =  0  (6) 

CD;  (1  +  QA)  +  kQA  =  0  (7) 
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(8) 


Substitution  of  o>,  from  Eq.  (7)  into  Eq.  (6)  yields 

cor2  =  We'1  k(l  -  k2)  b  +  k2a  =  0  , 


where 


a  =  QA/(1  +  QA)2  ,  b  =  B/(l  +  QA). 

The  extrema  of  cor  in  Eq.  (8)  must  satisfy 

We'1  (1  -  3k2)b  +  2ka  +  We'!k  (1  -  k2)  b'  +  kY  =  0.  (9) 

where  primes  denote  differentiation  with  k.  For  given  We  and  Q,  the  solution  of  (9)  gives  km 
the  wave  number  corresponding  to  the  most  amplified  disturbances.  Substitution  of  km  for 
k  in  Eq.  (8)  yields  (co^.  Substituting  this  value  for  lq  in  Eq.  (3),  we  have 

L  =  to/1  Cn  (a/ao)  (10) 


RESULTS 

The  dimensionless  intact  length  L  obtained  from  Eq.  (10)  with  fnCa/a^)  =  1  is  plotted  against 
the  dimensionless  jet  speed  U  =  WeQ  in  Fig.  1,  for  three  different  values  of  Q.  The  Rayleigh 
regime  I  where  QWe  <  0.8,  the  first  wind-induced  regime  II  where  0.8  <  QWe  <  1.2,  the 
second  wind-induced  regime  HI  where  1.2  <  QWe  <  5  and  the  atomization  regime  IV  are  also 
indicated  int  he  same  figure.  These  parameter  ranges  are  consistent  with  the  results  obtained 
in  the  author's  previous  investigations  [1-3],  on  the  mechanism  of  jet  breakup.  The  dotted 
lines  near  U  =  0  correspond  to  the  dripping  jet.  The  dripping  jet  arises  from  absolute 
instability  [8].  Some  experimental  results  taken  from  sterling  and  Sleicher  for  large  Reynolds 
number  are  also  given  in  the  same  figure  for  comparison.  It  should  be  pointed  out  that  the 
intact  length  discussed  here  refers  to  the  intact  liquid-gas  interface.  It  does  not  refer  to  the 
intact  core  length  [4,9],  the  determination  of  which  requires  nonlinear  theory. 

CONCLUSION 

Five  regimes  of  jet  breakup  have  been  recognized  from  the  solution  of  a  single  equation. 
These  regimes  are  dripping  jet,  the  Rayleigh  capillary  pinching,  first  wind-induced  breakup, 
second  wind-induced  breakup,  and  atomization.  The  dripping  jet  regime  occurs  when  the 
Weber  number  is  smaller  than  the  critical  Weber  number  [9].  The  Rayleigh  regime  occurs 
when  QWe  <  0.8.  The  first  wind-induced  regime  takes  place  in  0.8  <  QWe  <  1.2.  The 
second  wind-induced  regime  occurs  in  1.2  <  QWe  <  10.  Both  the  first  and  second  wind- 
induced  breakup  are  due  to  pressure  fluctuation  at  the  interface.  Finally,  the  atomization 
regime  starts  when  QWe  >  10.  The  present  classification  of  the  breakup  regimes  refers  to  the 
jet  interface,  but  does  not  refer  to  the  jet  core.  The  classification  of  the  regimes  of  jet  core 
breakup  requires  nonlinear  theories. 
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NOMENCLATURE 
a  ,  disturbance  amplitude 
k  ,  complex  wave  number 
L  ,  intact  length 
Q  ,  Density  ratio 
rD  ,  nozzle  radius 
Re  ,  Reynolds  number 
S  ,  Surface  tension 
We  ,  Weber  number 

Greek  symbols 

X  ,  argument  of  Bessel  function 
03  ,  complex  wave  frequency 
v  ,  kinematic  viscosity 

Subscripts 
0  ,  initial  value 

1  ,  liquid  phase 

2  ,  gas  phase 

i  ,  imaginary  part 
r  ,  real  part 
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Restoration  of  Laser 
Droplet  Slicing  Images 

Jingyi  Zhang  and  L.A.  Melton 
Department  of  Chemistry 
University  of  Texas  at  Dallas 
Richardson,  TX  75083-0688 

and 

Novel  Droplet  Slicing  Techniques 
Michael  Winter 

United  Technologies  Research  Center 
East  Hartford,  CT  06108 

In  droplet  slicing  imaging,  a  technique  in  which  the  midplane  of 
a  sub-millimeter  droplet  is  illuminated  with  a  laser  sheet  and 
the  fluorescence  is  viewed  (usually)  at  90°  to  the  sheet,  the 
image  is  significantly  distorted,  with  the  central  portion  being 
magnified  and  the  outer  portion  being  compressed  at  the  droplet 
edge.  At  UTD,  algorithms  and  programs  for  the  removal  of  these 
distortions  have  been  developed  and  tested. 

The  numerical  restoration  process  involves  several  steps. 

(1)  The  mapping  which  takes  the  actual  fluorescence  pattern 
within  the  droplet  (object  image)  to  that  measured  by  the  CCD 
camera  (camera  image)  is  constructed.  (2)  The  inverse  mapping, 
which  takes  the  camera  image  to  the  restored  object  image,  is 
constructed.  (3)  These  mappings  are  used  on  numerically 
constructed  images  to  test  the  degree  to  which  the  restored 
object  image  matches  the  object  image.  (4)  For  experimental 
data,  it  is  necessary  to  include  steps  involving  the  conversion 
from  Cartesian  coordinates  to  polar  coordinates  and  vice  versa. 

The  image  restoration  process  is  highly  unstable 
numerically.  While  it  is  possible  to  restore  simulated  90 
droplet  slicing  images  which  have  extremely  high  signal-to-noise 
ratios,  when  low  levels  of  random  noise  (1%)  or  perturbations  in 
the  index  of  refraction  are  added,  the  numerical  results  become 
meaningless  because  of  unconstrained  oscillations.  .The  numerical 
instability  can  be  suppressed  significantly  by  fitting  the  object 
image  to  a  set  of  basis  functions  and  performing  a  least  squares 
fit  for  the  best  coefficients  in  this  expansion. 

The  best  way  to  recover  the  image  all  the  way  to  the  edge  of 
the  droplet  is  to  carry  out  measurements  at  approximately  60  to 
the  laser  sheet.  In  this  manner,  the  image  on  the  acute  side  is 
recovered  easily  because  the  numerical  restoration  is  stable, 
while  the  information  on  the  oblique  side  is  lost  completely. 

AT  UTRC,  aerodynamic  droplet  generators  have  been  used. (1) 
to  reveal  the  interactions  between  internal  liquid  motions  in  the 
droplet  and  the  motions  in  the  surrounding  vapor  and  (2)  to 
generate  "layered  droplets"  in  which  fluorescent  acetone  vapor 
condenses  on  a  decane  droplet  and  reveals  the  streamlines  through 
its  fluorescence  (bright  on  a  dark  background) . 
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GEOMETRICAL  OPTICS  RAY  TRACING 


LEMS  IMAGE 

PLANE  PLANE 


COORDINATE  SYSTEMS 


COMPUTED  RADIAL  COORDINATES  MAPPING 
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MATHEMATICAL  FORMULATIONS 


IMAGE  DISTORTION  AND  RESTORATION 
90  DEGREE  OBSERVATION 


OBJECT  IMAGE  RESTORATION 


OFF  -  90  DEGREE  OBSERVATIONS 


NUMERICAL  SIMULATION  OF  OFF-90  DEGREE 

OBSERVATIONS 


IMAGE  DISTORTION  AND  RESTORATION 
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OBJECT  IMAGE  RESTORATION 


DROPLET  -  SLICING  IMAGE  AND  RESTORATION 


INTERNAL  STRUCTURE  OF  A 


The  Effect  of  Turbulence  on  Vaporization  and  Mixing  in  Fuel  Sprays 

D.  A.  Santavicca 

Department  of  Mechanical  Engineering 
Penn  State  University 
University  Park,  PA  16802 

The  objective  of  this  research  is  to  obtain  an  improved  understanding  of  the  behavior 
of  droplets  in  vaporizing  sprays,  particularly  under  conditions  typical  of  those  in  high  swirl, 
direct  injection  Diesel  engines.  Experiments  are  conducted  in  a  high  pressure,  high 
temperature  optically-accessible  flow  system  which  is  capable  of  operation  at  pressures  up  to 
70  atm,  temperatures  up  to  600  K,  gas  velocities  up  to  20  m/sec  and  turbulence  intensities  up 
to  40%.  Single  droplets,  50  to  500  micron  in  diameter,  are  produced  by  an  aerodynamic 
droplet  generator  and  transversely  injected  into  the  flow.  Measurements  are  made  of  the 
droplet  position,  size,  velocity  and  temperature  from  which  droplet  drag,  dispersion,  heating, 
vaporization  and  breakup  are  characterized. 

To  date,  the  main  results  from  this  study  are  the  following: 

1.  Under  laminar  flow  conditions,  vaporization  was  found  to  reduce  droplet 
drag,  in  quantitative  agreement  with  the  drag  correlation  of  Chiang  and 
Sirignano. 

2.  Under  laminar  flow  conditions,  droplet  drag  was  not  affected  by  unsteady 
curvilinear  motion. 

3.  Under  laminar  flow  conditions,  unsteady  curvilinear  motion  was  found  to 
result  in  small  but  non-negligible  droplet  lift  (CL/CD  ~  0.1),  but  only  at 
relatively  high  droplet  Reynolds  numbers  (20  <  Re  <  38). 

4.  Turbulence  was  found  to  result  in  an  apparent  increase  in  droplet  drag, 
however,  this  can  be  accounted  for  by  properly  defining  the  mean  relative 
velocity. 

5.  Turbulence  was  found  to  result  in  an  apparent  decrease  in  the  critical  Weber 
number  for  secondary  droplet  breakup,  however,  this  can  be  accounted  for  by 
properly  defining  the  mean  relative  velocity. 

6.  The  phenomenological  nature  of  secondary  breakup  was  observed  to  be 
fundamentally  different  in  turbulent  and  laminar  flows. 
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OBJECTIVE 


•  TO  CHARACTERIZE  THE  BEHAVIOR  OF  LIQUID  FUEL 

DROPLETS  IN  HIGHLY  CONVECTIVE  FLOWS  TYPICAL  OF 
THOSE  IN  DIESEL  FUEL  SPRAYS. 


•  TO  CHARACTERIZE  AND  QUANTIFY  DROPLET  DRAG, 

DISPERSION,  HEATING,  VAPORIZATION,  DISTORTION  AND 
BREAKUP. 


•  OF  PARTICULAR  INTEREST  ARE  THE  EFFECTS  OF 
TURBULENCE  AND  SUPERCRITICAL  CONDITIONS. 
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APPROACH 


•  EXPERIMENTS  ARE  CONDUCTED  WHERE  SINGLE  DROPLETS 
ARE  TRANSVERSELY  INJECTED  INTO  A  LAMINAR  OR 
TURBULENT,  SUB-  OR  SUPER-CRITICAL  FLOW. 


•  MEASUREMENTS  ARE  MADE  OF  THE  DROPLET  POSITION, 
TEMPERATURE,  SIZE,  SHAPE  AND  CONCENTRATION,  AS  A 
FUNCTION  OF  TIME,  FROM  WHICH  DROPLET  DRAG,  LIFT, 
DISPERSION,  HEATING,  VAPORIZATION,  DISTORTION  AND 
BREAKUP  ARE  DETERMINED. 
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HTGH  PRESSURE.  HIGH  TEMPERATURE 
TURBULENT  FLOW  SYSTEM 


Liqtid  Fuei 


M/ 


OPERATING  RANGE 


PRESSURE 
TEMPERATURE 
GAS  VELOCITY 
TURBULENCE  INTENSITY 
DROPLET  SIZE 


UP  TO  70  ATM 
UP  TO  600  K 
UP  TO  20  M/SEC 
UP  TO  40% 

50  TO  500  MICRONS 
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DROPLET  PARAMETER  MEASUREMENT  TECHNI 


STATUS 


•  RAMAN  IMAGING  OF  JET  BREAKUP  IN  A  SUPERCRITICAL 
ENVIRONMENT 

•  CALIBRATION  OF  DROPLET  THERMOMETRY 

•  VISUALIZATION  OF  DROPLET  VAPOR  WAKE 

RESULTS 

•  EFFECT  OF  VAPORIZATION  ON  DROPLET  DRAG  (LAMINAR) 

•  EFFECT  OF  UNSTEADY  CURVILINEAR  MOTION  ON  DROPLET 
DRAG  AND  LIFT  (LAMINAR) 

•  EFFECT  OF  TURBULENCE  ON  DROPLET  DRAG 

•  EFFECT  OF  TURBULENCE  ON  SECONDARY  BREAKUP 

•  EFFECT  OF  TURBULENCE  ON  DROPLET  DISPERSION 
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LAMINAR  EXPERIMENTS  - 


THE  EFFECT  OF 
VAPORIZATION  AND 
UNSTEADY 

CURVILINEAR  MOTION 


EXPERIMENTAL  CONDITIONS 


Exp.No. 

Dinit(^m) 

Tair  (K) 

P  (atm) 

BUBM 

Re 

®steadv 

1 

90 

293 

2.7 

1.25 

12-25 

0.2 

8 

90 

423 

4.4 

1.40 

10-25 

0.6 

11 

90 

473 

4.4 

180 

15-38 

0.9 

16 

103 

523 

5.76 

115 
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1.2 

DATA  ANALYSIS 


m- 


dc;^=-|p J0,-3M<P  l-CL(U3,y-Ud,y)  *CD(UgiX-(JafX) ) 


dt  2K»'  3  "  4 


m 


da^£.ifg\a-v^d2  lcLW,.x-ad.x)  *c„Wg.y-ad,y)  1  -mg 


~3t  2r’'  ’  4 


EFFECT  OF  VAPORIZATION  ON  DROPLET  DRAG 


EFFECT  OF  UNSTEADY  CimVHJNEAB  MOTION  ON  DROPLET  LIFT 
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DROPLET  TRAJECTORIES  IN  TURBULENT  FLOWS 


0.0000.0030.0050.008  0.0100.0130.0150.0180.020 


LOW  DISPERSION  CASE 


HIGH  DISPERSION  CASE 
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npnpi  FT  VFJ  nrTTY  IN  A  T , AMINAR  FLOW 


DROPLET  VFT  nriTY  IN  a  TURBULENT  FLOW 
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Computations  of  Sprays  in  a  Very  High-Pressure,  Constant  Volume  Chamber 

J.  Abraham 

Mechanical  Engineering  Department 
University  of  Minnesota 
Minneapolis,  MN  55455 

In  recent  years,  several  Diesel  engine  manufacturers  have  initiated  efforts  to  explore  ways  of 
achieving  increased  power-density  in  engines  through  the  use  of  extremely  high  cylinder 
pressures  (up  to  30  MPa).  The  engine  configurations  and  geometries  being  considered  are 
unique  and  it  is  apparent  that  rapid  progress  can  only  be  made  if  models  for  the  m-cylinder  flows, 
sprays  and  combustion  are  used  to  optimize  the  engine  geometry.  However,  the  accuracy  of  such 
models  have  not  been  assessed  at  the  conditions  being  proposed.  In  this  work,  comparisons  of 
computed  and  measured  penetrations  for  vaporizing  and  non-vaponzmg  sprays  m  a  fixed-volume 
chamber  will  be  presented.  The  ratio  of  ambient  gas  density  to  liquid  density  vanes  m  the  range 
0  005  -  0  25  It  will  be  shown  that  adequate  agreement  of  computed  and  measured  penetrations 
may  be  obtained  if  sufficiently  high  numerical  resolution  is  used  It  is  also  shown  that  as  the 
ratio  of  chamber  density/liquid  density  is  increased,  the  flammable  fraction  of  the  vapor  fuel 
increases. 
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OBJECTIVE: 


To  assess  the  accuracy  of  a  3-D  model  for  flows,  sprays  and  combustion  in  diesel  engines  under 
very  high-pressure  (density  of  ambient  gas/density  of  liquid  is  relatively  very  high)  engine 
operating  conditions. 


THE  MODEL: 

Numerics 

KIVA  (Los  Alamos) 

Physics 

1.  Turbulence  -  k-e 

2.  Boundary  layers  -  wall  functions 

3.  The  spray  model  -  (Princeton) 

(a)  atomization 

(b)  line  source 

(c)  turbulence  dispersion 

(d)  drop  breakup  (O'Rourke  &  Amsden,  Los  Alamos) 

(e)  collisions  &  coalescence 

4.  Autoignition  model 

5.  Combustion  model 
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EXPERIMENTAL  DATA:  (SANDIA  NATIONAL  LABORATORY) 


f  Non-Vaporizing  Sprays 

Case 

Pair/p  fuel 

TC(K) 

Pfuel  (MPa) 

A 

0.005 

300 

142 

B 

HHEXOEfll 

300 

145 

C 

0.033 

300 

142 

D 

0.047 

300 

139 

E 

0.061 

300 

142 

F 

300 

143 

G 

MPTHTTITP 

300 

143 

H 

0.103 

300 

139 

I 

0.148 

300 

142 

J 

0.243 

300 

138 

Vaporizing  Sprays  1 

Case 

Pair/pfuel 

TC(K) 

Pfuel  (MPa) 

A 

0.019 

900 

138 

B 

0.032 

899 

131 

C 

0.046 

908 

aa^pp— i i 

D 

0.061 

ns 

142 

E 

0.075 

mm 

145 
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PENETRATION  (MM) 


Measurements 


Measurements 


PENETRATION  VS  TIME:  VAPORIZING  SPRAYS 
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EFFECT  OF  GRID  RESOLUTION 
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NON-VAPORIZING  SPRAYS:  PENETRATION  VS  TIME 
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PENETRATION  (MM) 


NON-VAPORIZING  SPRAYS:  PENETRATION  VS  TIME 
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NON-VAPORIZING  SPRAYS:  PENETRATION  VS  TIME 
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PENETRATION  (MM) 


VAPORIZING  SPRAYS:  PENETRATION  VS  TIME 


TIME  (MS) 
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PENETRATION  (MM) 


CONCLUSIONS! 

1.  Comparisons  of  computed  and  measured  penetration  rates  of  sprays  injected  into  high- 
pressure  ambient  gas  show  adequate  agreement  provided  sufficient  numerical  resolution  is 

used. 

2.  As  chamber  pressure  density  is  increased,  the  flammable  fraction  of  the  vapor  fuel  increases 
and  the  rich  fraction  decreases.  If  this  persists  during  combustion,  the  combustion  will  be 
faster  as  the  chamber  density  is  increased. 
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Spray  Research  at  the  Engine  Research  Center 

P.  V.  Farrell 
with 

G.  L  Borman,  D.  E.  Foster,  J.  K.  Martin,  R.  D.  Reitz,  C.  J.  Rutland  and  lots  of  graduate 

students 

Engine  Research  Center 
University  of  Wisconsin-Madison 
Madison,  WI 53706 

A  large  number  of  two-phase  flow  and  related  spray  projects  are  underway  at  the 
University  of  Wisconsin-Madison  Engine  Research  Center.  For  diesel  engines,  the  spray 
projects  are  aimed  at  issues  which  include:  increasing  power  density  (power/package 
size);  examining  cold  start  spray  effects;  studies  of  fuel  effects  on  spray  behavior;  sprays 
in  some  special  application  engines;  and  studies  of  spray  effects  on  emissions.  For  spark 
ignited  engines  spray  work  emphasizes  spray  effects  on  mixing. 

At  the  Engine  Research  Center  (ERC)  a  variety  of  methods  are  employed  to 
investigate  the  issues  outlined.  These  methods  include  an  extensive  range  of 
experiments,  from  engine  experiments  to  specialized  off-engine  investigations  of  specific 
phenomena,  and  a  large  effort  in  applying  and  improving  computer  models  of  spray 
events.  The  diesel  injectors  studied  include  moderate  pressure  injectors  (500  bar),  high 
pressure  experimental  injectors  (1000  -  2000  bar),  and  high  pressure  commercial  or 
production  injectors  (1500  bar). 

The  experimental  tools  currently  applied  to  spray  studies  include  conventional 
measurements  (mass  flow  rate,  needle  lift,  pressures,...),  optical  methods  (photography, 
PDPA,  diffraction  based  sizing,  LIF,...),  and  some  developing  optical  methods  (PIV  and 
DSIV). 

The  modeling  tools  center  around  KIVA  as  the  core  of  the  model,  with  emphasis 
for  the  spray  on  the  submodels  most  connected  to  the  spray  predictions.  The  submodel 
developments  generally  parallel  experimental  efforts  in  the  same  area.  Modelers  can 
directly  tell  experimenters  what  kinds  of  data  they  need,  and  experimenters  can  discuss 
with  modelers  what  can  be  delivered.  A  very  important  part  of  the  ERC  activities 
involves  this  interaction  between  model  development  and  validation  and  experimental 
focus,  design,  and  operation.  We  think  it  is  vital  that  students  and  faculty  recognize  the 
valued  requirements,  and  limitations  of  each  of  these  approaches  in  order  to  progress  on 
all  fronts. 

This  presentation  will  attempt  to  outline  the  range  of  projects  currently  underway 
at  the  ERC  and  to  indicate  the  general  objectives  of  the  projects  and  what  we  hope  will  be 
the  impact  of  these  projects.  Many  of  the  projects  described  are  funded  by  ARO, 
however  for  completeness,  some  non-ARO  funded  projects  will  also  be  mentioned.  In 
addition  to  the  project  survey,  two  areas  of  projects  will  be  examined  in  somewhat  more 
detail  in  order  to  present  a  little  more  detail  for  these  topics  and  some  results. 
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Methods  for  Spray 

Studies  at  ERC  Methods  for  Spray 
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Particulate  (g/bhp-hr) 


Fuel  System 

Injector  Type  Common  Rail  -  Electronically  Controlled 

independent  control  of  the  duration  and  timing  of 
up  to  4  injections  per  cycle 

Injection  Pressure  (variable)  90  MPa  (13,230  psia) 

Number  of  nozzle  orifices  6 
Nozzle  diameter  0.26  mm 

Spray  Angle  (included)  125' 

Nozzle  Cd  0.78 

Fuel  Type  Amoco  1994  Low  Sulfur  Emissions  Certification  Fuel 


m  Pulst:  Degrees 

I  I 

7-(2)-44-{10)-49 


Particulate  vs.  NOx  75%  Load  Particulate  vs.  NOx  75%  Load 


NOx  (g/bhp-hr)  NOx  (g/bhp-hr) 


BSFC  and  Particulate  at  NOx  =  5  g/bhp-hr 
75%  Load  Cases 


[Model  based  on  Jin  & 
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Motor 
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ITEX  IMAGING  BOARD 
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Table  4.1.  Injector  Calibration  Parameters 


Parameter 

Lombardini 

specifications 

Large  nozzle 
D-2  fuel 

Large  Nozzle 
Jet-A  fuel 

Small  Nozzle 
D-2  Fuel 

Injected  mass 
(mg) 

=20 

19.7 

19.7 

7.7 

Nozzle  Open. 
Press.  (Mpa) 

19-21 

19.0 

19.0 

19.0 

Mean  Injection 
Press.  (Mpa) 

25.1 

25.4 

25.1 

Peak  Injection 
Press.  (Mpa) 

33.8 

33.8 

40.4 

Injection 
Duration  (ms) . 

1.7 

1.7 

1.8 

Mean  Inj. 
Velocity  (m/s) 

186.3 

192.6 

198.7 

Table  4.3.  System  Specifications  for  the  Laser  Interrupt  Test  and 
Variables  Studied  with  Ranges. 


Injection  System  Specification: 

Lucas  CAV  injector  with  Bosch  injection 
PumP- _ _ _ 

Pump  Speed 

1500  RPM 

Radial  Penetration  Distances  Rf(mm) 

10,  15,  20,  30 

Nozzle  Diameter,  d  (mm) 

0.180  (+),  0.306  (-) 

Impingement  Distance,  L,  (mm) 

6,(0  38  (+) 

Target  Diameter,  Df  (mm) 

14.7  (0,18.3  (+) 

Target  Temperature,  T  (°C) 

150  (-) ,  315  (+) 

Ambient  Density,  pa  (kg/m3) 

■14.34(0,  22.07  (+) 

Ambient  Pressure,  Pa  (Mpa) 

1.23  /  1.93  No  (-).  1.51/2.31  80%  N?(+) 

Target  Materials  Conductivity  • 

(W/m2K  (5)  400K) 

16.6(0,197  (+) 

Fuel  Type 

Amoco  Low  Sulfur  D-2,  Commercial  Jet- 
A 
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pray  image  film  frame  6.  (b)  and  (c)  Super-imposed  edged  spray 
spray  and  cropped  spray  development  diagrams. 


Converging  Lenses 


FWIrf  PwmJm  (mm) 


Tm  (rm)  Tmw  (m»)  ("»)_ 


Figure  6.40.  (a-c)Comparisons  of  spray  penetration  of  a  normal  cone  spray  and 
a  spray  impinged  against  a  raised  target 


Jtm  (m») 


Figure  6.38.(a-c).  Comparisons  of  penetration  of  a  spray  impinged  upon  a  flat 
wall  and  a  spray  impinged  against  a  raised  target 
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Drop  Impinges 
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The  predicted  spray  pattern  for  comparison  to  Senda  et  al. 
(1991)  at  1  ms.  2  ms.  3  ms,  6  ms,  and  9  ms  after  injection  (T  =  500  K) 


Advanced  Diesel  Injection  Strategies 

F.V.  Bracco 
Princeton  University 

Direct  injection  of  natural  gas  in  Diesel  engines  by  several  engine  companies  has  yielded  the 
surprising  result  that  combustion  time  is  much  longer  with  natural  gas  than  with  Diesel  fuel  (the 
ignition  delay  time  is  not  much  longer).  Natural  gas  also  yields  unexpectedly  high  soot  emission. 
These  trends  are  surprising  because  it  has  always  be  assumed  that  the  faster  the  vaporization  the 
faster  the  combustion  and  the  limit  of  very  fast  vaporization  is  the  gas  jet 

Results  are  presented  of  3-D  computations  of  direct  injection  of  gaseous  methane  and  of 
liquid  tetradecane  through  a  multi-hole  injector  into  a  Diesel  engine.  The  study  focusses  on  the 
distribution  of  fuel/air  ratio  within  the  resulting  gas  and  spray  jets  under  typical  Diesel  conditions 
prior  to  ignition.  It  is  shown  that  for  a  significant  time  after  start  of  injection,  the  fraction  of  the 
vapor  fuel  which  is  in  richer-than-flammable  mixtures  is  greater  in  gas  jets  than  in  sprays.  For 
methane  injection,  it  is  also  shown  that  changing  some  of  the  flow  conditions  in  the  engine  or 
going  to  a  poppet-type  injector,  does  not  result  in  improved  mixing.  It  is  concluded  that  in  general 
a  spray  mixes  faster  than  a  gas  jet  in  direct-injection  Diesel  engines  and  possible  explanations  for 
this  unexpected  result  are  provided.  An  experiment  is  also  in  progress  to  check  the  proposed 
explanations. 
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Fig.2:  Velocity  flowfield  and  equivalence  ratio  contour 
plots  for  the  baseline  case. 
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Fig.5:  Velocity  flowfield  and  equivalence  ratio  contour 
plots  for  the  case  with  swirl. 
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Fig. 6:  ”f  -  6n  plot  for  the  case  with  swirl. 
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Fig. 7:  Velocity  flowfield  and  equivalence  ratio  contour 
plots  for  the  case  with  a  15°  angle  between  the  engine 
head  and  injector  hole  orientation. 


Fig.8:  ”f  -  6"  plot  for  the  case  with  a  15°  angle  be¬ 
tween  the  engine  head  and  injector  hole  orientation. 
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Fig.9:  Velocity  flowfield  and  equivalence  ratio  contour 
plots  for  the  case  with  increased  squish. 


EQUIVALENCE  RATIO 

A 

Fig.10:  ”f  -  6 ”  plot  for  the  case  with  increased  squish. 

Fig.ll:  Velocity  flowfield  and  equivalence  ratio  con¬ 
tour  plots  for  gas  injection  into  a  cylindrical  bowl  with  a 
conical  gas  jet. 


15-6 


Fig.  12:  "f  -  6"  plot  for  gas  injection  into  a  cylindrical 
»wl  with  a  conical  gas  jet  (reference  case). 


:o 


Fig.  13:  ”f  -  0”  plot  for  gas  injection  into  a  cylindrical 
>wl  with  a  conical  gas  jet  and  with  the  turbulence  in- 
nsity  increased  by  a  factor  of  2  over  the  reference  case. 


o 


Fig.  14:  ”f  -  9V  plot  for  gas  injection  into  a  cylindrical 
bowl  with  a  conical  gas  jet  and  with  the  turbulence  length 
scale  increased  by  a  factor  of  2  over  the  reference  case. 


Fig.  15:  Velocity  flowfield  and  equivalence  ratio  con¬ 
tour  plots  for  gas  injection  into  a  fixed  volume,  2ms  after 
the  start  of  injection. 


Fig.16:  Velocity  flowfield  and  equivalence  ratio  con¬ 
tour  plots  for  direct  injection  of  liquid  into  a  fixed  volume 
2ms  after  the  start  of  injection. 


volume. 


the  fixed  volume  with  the  Sauter  Mean  Radius  of  the 
initial  drops  increased  by  a  factor  of  10. 


Fig.21:  Fraction  of  injected  fuel  momentum  in  fuel 
(curve  A),  ambient  gas  (curve  B)  and  in  the  pressure 
difference  (curve  C)  at  3.2cm  from  the  injector  for  gas 
injection  into  the  fixed  volume. 


the  fixed  volume. 


initial  drops  increased  by  a  factor  of  25. 


Fig.22:  Fraction  of  injected  fuel  momentum  in  fuel 
(curve  A),  ambient  gas  (curve  B)  and  in  the  liquid  phase 
(curve  C)  at  3.2cm  from  the  injector  for  liquid. injection 
into  the  fixed  volume. 
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Fig. 23:  Penetration  of  gas  jets  for  gas  injection  into  Fig.24:  ”f  -  9"  plot  for  gas  jet  in  an  infinite  environ- 


an  infinite  environment. 
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High-Speed  Four-Color  Infrared  Digital  Imaging  for  Studying 
In-Cylinder  Processes  in  a  DI  Diesel  Engine 


K.  T.  Rhee 

Department  of  Mechanical  and  Aerospace  Engineering 
Rutgers,  The  State  University  of  New  Jersey 
Piscataway,  NJ  08855 


Four  geometrically  identical  digital  IR  images  of  in-cylinder  events  are  simultaneously  obtained 
in  respective  spectral  bands,  which  is  individually  triggered  by  the  crank-angle  signal  generated 
from  the  host  direct  injection  (DI)  diesel  engine. 

A  new  spectrometric  information  processing  algorithm  is  developed  in  order  to 
determine  distributions  of  temperature  and  species  in  the  combustion  chamber  using  the 
measurements. 

Apparatus  Construction.  The  improved  engine  apparatus  having  optical  access 
has  recently  been  incorporated  with  a  new  electronic  controlled  fuel  injection  system  fabricated 
in  our  laboratory.  This  system  at  present  delivers  a  injection  pressure  of  as  high  as  170  MPa 
(25,000psi),  which  will  be  increased  to  well  over  200MPa  in  the  next  system. 

While  our  early  two-color  imaging  system  was  employed  in  engine  studies, 
development  of  the  four-color  spectral  imaging  system  has  been  performed  at  the  same  time. 

This  includes:  (a)  a  newly  developed  reflective  optical  arrangement  having  three 
custom-made  spectral  beam  splitters;  (b)  four  IR  imaging  camera  heads;  (c)  new  electronic 
packages  having  printed  boards  (unlike  the  hand-wired  electronic  boards  employed  in  the 
previous  systems)  for  driving  the  camera  heads  and  simultaneously  retrieving  the  output  from 
the  cryogenically  cooled  imagers;  and  (d)  several  new  computer  programs  for  image  processing 
and  display. 

Experimental  Results.  The  present  ARO-sponsored  study  has  been  joined  by  Ford 
Motor  Co.  and  Texaco  Research  Center  in  order  to  study  SI  engine  combustion  by  using  this 

new  diagnostic  system.  .  . 

(1)  SI  engine  combustion  has  been  the  subject  to  initially  study  by  using  this  new  tool 
because  its  in-cylinder  environment  is  easier  to  manage  than  its  counter  part,  the  Cl  engine. 

This  co-sponsored  work  is  being  directed  to  investigation  of  several  engine  events,  including 
the  knock  processes,  flame  kernel  growth,  deposit  image,  and  cold  start. 

(2)  The  most  significant  finding  from  the  DI  diesel  engine  study  reported  in  the 
previous  period  was  the  low  temperature  inferred  in  the  regions  near  the  axis  of  the  spray 
plume,  which  has  been  further  studied  in  order  to  explain  the  reasons. 

Among  the  possible  reasons  are:  The  fuel/air  in  the  regions  are  at  speeds  too  high  to 
permit  stable  formation  of  reaction  fronts,  which  seems  to  be  supported  by  our  video  movie 
produced  using  the  new  electronic  fuel  injector  and  the  four-color  system.  The  other  reason  is 
that  the  fuel  vaporization  results  in  lower  temperatures  than  required  for  self-ignition  in  the 
regions.  The  next  possible  reason  is  that  the  combination  of  the  above  effects.  That  is,  even 
though  the  mixtures  attain  either  marginally  or  sufficiently  high  temperatures  for  self-ignition,  it 
will  be  difficult  to  produce  established  reaction  fronts  when  the  reactive  mixture  is  at  high 
speeds  along  the  axis. 

Future  Activities.  Since  our  imaging  system  is  now  operational,  which  we  will 
continue  to  improve,  we  will  obtain  measurements  and  analyze  them  using  our  new  process 
software  in  order  to  implement  the  project  goal. 

The  variables  to  be  immediately  studied  include  different  injector  tips  provided  by 
Cummins  Engine  and  new  fuels  expected  from  the  DOE  Alternative  Fuels  Utilization  Program. 
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o.  Objectives 
o.  Methods 

o.  Analysis  of  Methods 
o.  Progress 
o.  Future  Activities 


1.  Develop  a  New  High-Speed  Multispectra  IR  Imaging  System 

for  Engine  Studies. 

2.  Develop  Spectrometric  Data  Analysis  Algorithm. 

3.  Investigate  In-Cylinder  Processes  of  a  DI  Diesel  Engine 

during  the  Entire  Cycle  Period: 

.  Residual  Gas  Behaviors, 

.  Temporal  and  Spatial  Temperature  Distributions, 

.  Combustion  Product  Distributions, 

.  Heat  Transfer,  and  Others. 


W®(r& 

1.  A  New  Multispectra  Imaging  System. 

2.  Conventional  Two-Color  Method  for  Solid  Wall. 

3.  New  Dual-Band  Method  for  Gaseous  Mixtures. 

4.  A  New  Electronic-Controlled  Fuel  Injection  System. 

5.  Bench  Test  and  Engine  Experiment. 
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Fig.  1.  Assembly  Drawing  of  Four-Color  Imaging  System  with  Reflective  Optics  Arrangement 
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Fig.  2.  Flow-Chart  of  the  Dual-Band  Quantitative  Imaging. 
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Fig.  5.  Schematics  of  SI  Engine  Apparatus. 

(Co-Sponsors:  Ford  Motor  and  Texaco) 


Fig.  6.  Diesel  Engine  with  Optical  Access 
(ARO  Combustion  Research) 


Fig.  7.  Four  Color-Imaging  System  Lined  up  with  an  Engine  Apparatus. 


Temperature  Water  Vapor 


Fig.  8.  Instantaneous  Quantitative  Flame  Images  of  Hydrogen-Air  Burner  (Fig.  4). 
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Fig.  9.  High-Speed  IR  Images  of  Flame  Kernel  Growth  in  an  SI  Engine  (Fig.  5). 
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Fig  10.  High-Speed  Spectral  IR  Images  of  a  Cl  Engine  Combustion  with 

an  Electronic  Controlled  Fuel  Injector  (170  MPa,  Injection  Pressure). 
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Combustion  Instability  During  Cold  Starting  of  Diesel  Engines 

N.  A.  Henein,  M.  Lai 
Center  for  Automotive  Research 
Wayne  State  University 
Detroit,  MI  48202 

Combustion  instability  during  cold  starting  of  diesel  engines  causes  many  problems  such  as 
long  periods  of  cranking  and  the  emission  of  excessive  amounts  of  unbumed  fuel  in  the 
form  of  white  smoke,  stalling  after  an  instable  combustion  or  a  complete  failure  of  the 
engine  to  start.  Such  problems  are  of  major  concern  to  engine  manufacturers  and  specially 
to  the  Army  under  the  critical  field  operating  conditions. 

The  objectives  of  this  research  program  are:  i-  to  investigate  the  nature  of 
combustion  instability,  its  randomness,  and  dependence  on  the  fuel  properties  and  engine 
design,  ii-  to  investigate  the  causes  of  combustion  instability  and  iii-  to  propose  solutions. 

Tests  are  conducted  on  three  single  cylinder,  direct  injection,  four-stroke-cycle,  diesel 
engines  of  different  designs.  The  first  is  a  TACOM  research  engine,  water  cooled, 
connected  to  an  electric  dynamometer  for  motoring  and  power  absorption,  equipped  with 
an  endoscope.  The  second  is  an  industrial  stand-alone  engine  equipped  with  its  electric 
starter,  battery  and  fuel  tank.  Both  these  engines  are  installed  in  a  cold  room  where  the 
ambient  temperature  can  be  controlled  between  the  normal  room  temperature  and  -50  C. 
The  cold  starting  tests  are  conducted  after  the  engines  are  soaked  for  6  to  8  hours.  The 
third  engine  is  an  AVL  520  water  cooled  engine  modified  for  optical  access  through  the 
piston  top,  all  around  the  top  of  the  cylinder  and  through  the  cylinder  head.  The  engine 
is  connected  to  an  electric  dynamometer  for  motoring.  The  intake  air  temperature  is 
controlled.  Fuel  injection  occurs  only  for  two  consecutive  cycles  while  the  engine  is  being 
motored.  Traces  are  taken  for  the  fuel  pressure  before  the  injector  and  for  the  cylinder  gas 
pressure.  The  diagnostic  instrumentation  includes  the  following:  (a)  a  high  speed  movie 
camera  to  study  the  spray  parameters  such  as  penetration,  angle,  impingement  on  and 
reflection  off  the  walls  of  the  combustion  chamber,  start  of  the  ignition  process,  flame 
propagation  and  smoke  intensity,  (b)  Exciplex  LIF  to  study  spray  evaporation  (without 
combustion)  in  a  nitrogen  atmosphere,  (c)  High  speed  video  for  2-color  pyrometry  for 
simultaneous  flame  temperature  measurements,  and  OH  chemiluminescence  imaging,  (d) 
Lase  Incandescence  Imaging  (LII)  to  investigate  the  spatial  distribution  of  soot  in  the  spray. 

The  experimental  results  on  the  three  engines  indicated  the  following:  (a) 
Combustion  instability  is  not  a  random  process  and  combustion  may  occur  every  sixteen, 
twelve,  eight  or  four  strokes  depending  on  the  ambient  temperature,  (b)  Combustion 
instability  is  not  engine  or  fuel  specific,  (c)  Combustion  instability  is  caused  by  a 
combination  of  factors  related  to  engine  dynamics  and  combustion  kinetics,  during  the 
transient  cold  starting  process. 

High  speed  photographs  of  the  two  consecutive  cycles  were  obtained  at  different 
intake  air  temperatures:  23°C,  18°C  and  13°C.  At  23°C  the  engine  fired  and  combustion  was 
stable  in  the  two  consecutive  cycles.  At  18°C  the  engine  fired  in  the  first  cycle  and  misfired 
in  the  second  cycle.  At  13°C  the  engine  misfired  in  the  two  cycles. 

Future  plans  are  for  the  use  of  different  laser  based  techniques  to  investigate  the 

cause(s)  of  combustion  instability. 


17-1 


OVERALL  GOAL 


IMPROVE  THE  COLD  STARTING  OF  MILITARY 
DIESEL  ENGINES. 

BACKGROUND 

1.  DIESEL  COLD  STARTING  PROBLEMS  INCLUDE: 

.  LONG  PERIOD  OF  CRANKING  AND  THE  EMISSION 
OF  EXCESSIVE  AMOUNTS  OF  UNBURNED  FUEL, 
(WHITE  SMOKE). 

.  STALLING  AFTER  AN  INSTABLE  COMBUSTION. 

.  COMPLETE  FAILURE  OF  THE  ENGINE  TO  START. 

2.  MOST  OF  THE  PROBLEMS  ARE  RELATED  TO  AN 
INSTABLE  COMBUSTION  PROCESS. 
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OBJECTIVES 


TO  ANSWER  THE  FOLLOWING  QUESTIONS: 

1.  IS  COMBUSTION  INSTABILITY  DURING  COLD 
STARTING  A  RANDOM  PROCESS  ? 

2.  IS  COMBUSTION  INSTABILITY  CAUSED  BY  A 
FAILURE  OF  THE  AUTOIGNITION  PROCESS  ? 

3.  IS  COMBUSTION  INSTABILITY  SPECIFIC  TO  ONLY 
SOME  FUELS  ? 

4.  IS  COMBUSTION  INSTABILITY  SPECIFIC  TO  ONLY 
SOME  ENGINES  ? 

5.  WHAT  IS  THE  EFFECT  OF  TEMPERATURE  ON 
COMBUSTION  INSTABILITY  ? 

6.  WHAT  ARE  THE  CAUSES  OF  COMBUSTION 
INSTABILITY  ? 
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APPROACH 


THEORETICAL 

.  ANALYZE  THE  COMBUSTION  PROCESS  UNDER  THE 
TRANSIENT  THERMAL  AND  DYNAMIC  ENGINE 
CONDITIONS  DURING  COLD  CRANKING  AND 
STARTING. 

.  DEVELOP  A  MATHEMATICAL  MODEL  FOR  THE 
ENERGY  RELEASE  IN  SUCCESSIVE  ENGINE  CYCLES, 
CONSIDERING  THE  FOLLOWING: 

.  INSTANTANEOUS  ANGULAR  VELOCITY 

.  INSTANTANEOUS  BLOWBY 

.  INSTANTANEOUS  FUEL  INJECTION 

.  INSTANTANEOUS  BURNING  RATES 
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F.XPF.RTMENTAL 

F.NGTNES:  DIRECT  INJECTION,  4-STROKE-CYCLE 

(i)  TACOM,  DYNAMOMETER  DRIVEN,  WATER 
COOLED,  WITH  AN  ENDOSCOPE  FOR  IN 
CYLINDER  IMAGING. 

(ii)  DEUTZ,  STAND-ALONE,  AIR  COOLED. 

(iii)  AVL,  MODIFIED  FOR  OPTICAL  ACCESS, 
DYNAMOMETER  DRIVEN,  WATER  COOLED. 

.  HIGH  SPEED  MOVIE  CAMERA 

.  EXCIPLEX  LIF 

.  HIGH  SPEED  2-COLOR  PYROMETRY 
.  LASER  INCANDESCENCE  IMAGING 
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DETAILED  (P  -  0)  TRACES  FOR  THE  12-STROKE- 

CYCLEOPERATION 
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(P  -  6)  AND  (m  -0)  TRACES  FOR  DEUTZ  ENGINE, 
OPERATING  ON  8-STROKE-CYCLE 
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Fig  1 4:  Detailed  (P  -  9  )Traces  for  the  8-Stroke-cycle  Operation 


DETAILED  (P  -  0)  TRACES  FOR  THE  8-STROKE-CYCLE 

OPERATION 
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CYLINDER  GAS  TEMPERATURE  AT  START  OF  INJECTION 
FOR  CYCLES  22  TO  29  [DEUTZ  ENGINE] 


rad/ sec 
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COMBUSTION  INSTABILITY  DURING  COLD 
STARTING  IN  TACOM  ENGINE 


Bate  of  Heat  Release  |mg/CA*]  Rate  of  Heat  Re|ease  (mg/cA*] 


Cycle  23  To  23  Increment  1  Smoothing  Factor-0  [R.S.] 
Fuel  refl  TamD  0*C  Type:  STARING  FIRST  Timing:  23‘BTDG 


Crankshaft  Position  [°] 


Cyole  24  To  24  Increment  1  Smoothing  Fador-0  [R.S.J 
Fuel  reTl  TamD  Q«C  Type:  STAR l  NG  Fi  RST  Timing:  23*BTDC 


336  344  352  360  368  376  3i 

Crankshaft  Position  t°] 


RATE  OF  HEAT  RELEASE  FOR  CYCLES  23  AND  24- 


4 


MODIFIED  AVL  ENGINE:  EXTENDED 
CYLINDER/PISTON  AND  HEAD  ASSEMBLY 
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Pressure  Pressure  (bar)  chamber  pressure  (bar) 


EFFECT  OF  INLET  AIR  TEMPERATURE  ON 
COMBUSTION  STABILITY  [OPTIC  ACCESS  ENGINE] 
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CONCLUSIONS 


1-  COMBUSTION  INSTABILITY  DURING  COLD 
STARTING  IS  NOT  A  RANDOM  PROCESS.  ENGINE 
MAY  OPERATE  ON  A  16,  12,  8  OR  4-STROKE-CYCLE, 
DEPENDING  ON  THE  AIR  TEMPERATURE. 

2-  COMBUSTION  INSTABILITY  DURING  COLD 
STARTING  MAY  NOT  BE  CAUSED  BY  A  FAILURE  OF 
THE  AUTOIGNITION  PROCESS. 

3-  COMBUSTION  INSTABILITY  DURING  COLD 
STARTING  IS  NOT  ENGINE  OR  FUEL  SPECIFIC. 

4-  COMBUSTION  INSTABILITY  DURING  COLD 
STARTING  IS  CAUSED  BY  A  COMBINATION  OF 
FACTORS  RELATED  TO  ENGINE  DYNAMICS 
AND  COMBUSTION  KINETICS. 
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Effect  of  Thin  Ceramic  Coatings  on  Combustion  and  Emissions  in  a  DI  Diesel  Engine 


D.  E.  Klett 

Mechanical  Engineering  Department 
North  Carolina  A&T  State  University 
Greensboro,  NC  27411 

E.  M.  Afify 

Mechanical  and  Aerospace  Engineering  Department 
North  Carolina  State  University 
Raleigh,  NC  27607 


An  experimental  and  analytical  study  on  the  effect  of  thin  ceramic  coatings  on  soot  and  NOx 
emissions  and  performance  of  a  direct  injection  diesel  engine  is  being  conducted  jointly  between 
North  Carolina  A&T  State  University  and  North  Carolina  State  University. 

Performance  and  emissions  data  have  been  gathered  on  a  normally  aspirated  Ricardo 
Hydra  single  cylinder  DI  engine  with  various  combinations  of  ceramic  coatings  installed.  Thin 
ceramic  thermal  barrier  coatings  were  applied  to  the  piston  crown  and  bowl,  the  head  and 
valves,  and  the  cylinder  liner.  The  coated  piston  and  head  were  run  singly  and  in  combination 
with  the  cylinder  liner  to  investigate  the  effects  of  these  different  coated  surfaces  on  emissions 
and  performance  for  two  different  pure  hydrocarbon  fuels,  hexadecane  and  dodecane.  Coating 
the  piston  crown  alone  results  in  generally  lower  cylinder  pressure,  lower  brake  specific  fuel 
consumption  and  lower  NOx  emission  compared  to  the  baseline  engine.  Soot  emission  is 
typically  increased  below  2000  RPM  and  decreased  above  2000  RPM.  Coating  the  head  alone 
reduces  cylinder  pressure,  but  generally  increases  specific  fuel  consumption  and  NOx  and  soot 
emissions. 

The  analytical  portion  of  the  study  involved  modifications  to  the  KIVA-II  code  and  its 
use  to  model  the  Hydra  engine  with  the  thermal  coatings.  Modifications  to  the  code  include 
incorporation  of  an  eddy-break-up  combustion  model  to  replace  the  standard  Arrhenius  single 
reaction  model.  A  time  dependent  combustion  surface  temperature  was  also  incorporated  to 
simulate  the  effects  of  thermal  barrier  coatings  on  cylinder  temperature,  pressure  and  NOx 
production.  A  soot  model  was  added  to  the  code  following  the  work  of  Magnussen  and 
Hjertager.  The  EBU  model  gives  better  results  for  the  diffusional  portion  of  the  combustion 
process,  but  fails  to  adequately  model  the  premixed  combustion,  typically  resulting  in  a  lower 
predicted  peak  cylinder  pressure  (and  temperature)  than  predicted  by  the  single  reaction  model 
and  shown  by  experiments.  Consequently,  the  EBU  model  also  under-predicts  NO  emission  to 
a  greater  degree  than  the  SR  model.  The  KIVA-II  modeling  has  led  to  an  understanding  of  the 
effect  of  coating  the  piston  on  NO  production.  The  hotter  piston  crown  warms  the  intake  air, 
shortening  ignition  delay  and  decreasing  the  ratio  of  premixed  to  diffusion  combustion, 
ultimately  resulting  in  lower  peak  cylinder  temperature  and  reduced  NO.  The  KIVA-II  results 
agree  reasonably  well  with  the  experimental  data  for  cylinder  pressure  and  NO  and  soot 
emission. 

The  work  is  continuing  with  the  use  of  JP-8  fuel.  Performance  and  emissions  will  be 
compared  between  the  baseline  engine  and  with  the  various  coated  surfaces  with  JP-8  as  the  fuel. 
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•  INVESTIGATED  EFFECTS  OF  THIN  THERMAL  BARRIER 
COATINGS  ON  VARIOUS  SURFACES  ON  PERFORMANCE 
AND  EMISSIONS  OF  THE  RICARDO  HYDRA  DI  ENGINE 


•  TESTS  WERE  MADE  WITH  VARIOUS  COMBINATIONS  OF 
SURFACE  COATINGS 

•  PISTON  CROWN 

•  CYLINDER  HEAD  AND  VALVE  FACES 

•  PISTON,  HEAD  AND  CYLINDER  LINER 


•  EFFECT  ON  BSFC,  IGNITION  DELAY,  SMOKE,  NOX 

•  EXPERIMENTAL  RESULTS  COMPARED  WITH  KTVA-II 
MODELING  RESULTS 


•  KTVA-II  MODEL  INCLUDES  EDDY-BREAK-UP  COMBUSTION 
MODEL,  SOOT  MODEL  AND  COATED  SURFACE  THERMAL 
MODEL 


•  WORK  IS  CONTINUING  CURRENTLY  WITH  JP-8  FUEL 


THERMAL  BARRIER  COATINGS 


Piston  Crown  -  0.25  mm  slurry  sprayed  PSZ 

85%  Cubic  Zirconia,  calcium  partial  stabilization 
10%  Tungsten  Cobalt  Chrome  powder 
5%  Cr^Og 

Cylinder  Liner  -  0.75  mm  Plasma  Sprayed  YSZ  above  TBR 
0.2  mm  slurry  sprayed  PSZ  wear  coat  over 
entire  liner  surface 

Head  and  Valves  -  0.5  mm  TBC  with  hollow  Alumina  spheres 
65%  Silica 
15%  PSZ 

7%  Tungsten  Cobalt  Chrome  powder 

8%  CTgOg 

5%  Alumina  spheres  _ 


|  TEST  MATRIX 

lnj  Timing  (Deg  BTDC) 

16 

18 

20 

22 

BMEP  (BAR) 

2.24 

3.35 

3.91 

4.47 

Speed  (RPM) 

1000 

1500 

2000 

2500 

Engine  Build 

Baseline 

Coated 

Piston 

Coated 

Head 

Coated 

Piston, 

Head,  Liner 
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Af  138  DBTDC 


At  TDC 


Conipiilotionnl  Mejlt  Used  to  Model  Kienrdo  Hydro  Dl  Engine 
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1000  RPM  &  20  DBTDC 


1500  RPM  k  20  DB1DC 


300  <00  500  600 

BMEP  (KPo) 


400 

BMEP  (KPo) 


2000  RPM  &  20  DBTDC 


2500  RPM  k  20  DBTDC 


400 

BMEP  (KPo) 


400 

BMEP  (KPo) 


BSFC  Versus  Load  for  Dodecane  Fuel  at  20  Degrees  BTDC 


(iudd)  xqn  (ujdd)  XQN 


Soot  Concentration  (fl/m~3) 


1000  RPM  k  20  DBTDC 


1500  RPM  k  20  DBTDC 


(a) 


(b) 


2000  RPM  k  20  DBTDC 


2500  RPM  k  20  DBTDC 


(c) 


(d) 


Soot  Concentration  Versus  Load  for  Dodecane  at  20  DBTDC 
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NO  Emission  Versus  Cr.nk  Angle  For  Hes.dec.ne  Fuel 
1500  RFM,  16  Deg  BTDC,  213  KP»  BMEP 


NO  Emission  Versus  Cr.nk  Angle  Tor  l!e*«dec«ne  Fuel 
1500  RPM,  18  Deg  BTDC,  391  KP*  BMEP 


NO  Emission  Versus  Cr.nk  Angle  For  Dodec.ne  Fuel 
1500  RPM,  16  Deg  BTDC,  223  KP.  BMEP 


NO  Emission  Versus  Cr.nk  Angle  for  Dodec.ne  Fuel 
1500  RPM,  18  Deg  BTDC,  223  KP.  BMEP 
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SOOT  CONCENTRATION  (G/M**3) 


SOOT  EMISSION  (G/M**3)  S00T  EM|SS10N  (G/M**3) 


0.75 


0.5 


INJ.  TIMING=1B  DEG  BTDC  - 

FUEL:  HEXADECANE 

1500  RPM 

STP 


KIVA 

experiment 


0.25 


0 


1 


x 

3 


x 

4 


BMEP  (BAR) 

Soot  Emission  Versus  Load  Tor  Hexadecane 


Soot  Emission  Versus  Load  for  Dodecane 


0 


£ 

3 


550 

525 

500 

475 

450 

425 

400 

-360  -180  0  180  360 

Crank  Angle  (DEG) 

Assumed  Coaled  Surface  Temperature  Versus  Crank  Angle  Profile 


_  Hon— Ineirloted  Surfoce 

.  Insulated  Surface 

- 

_ 1 _ 1 _ 1 _ 

Experimental  Results  For  ISO  Emission  Versus  Load 
For  Hexadecane  at  1500  RFM  and  22  DBTDC  Injection  Timing 


BMEP  (KPa) 

KJVA  Results  For  NO  Emission  Versus  Load 
For  Ilexadecnne  at  1500  RTM  and  22  DBTDC  Injection  Timing 
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Cylinder  Gns  Pressure  Versus  Crank  Angle 


Cylinder  Gas  Temperature  Versus  Crank  Angle 
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Mass  of  Fuel  Burnt  Versus  Crank  Angle 


...,,,.7  pitton  buKArttdA 

_ _ food  kwvteUd  ♦ 
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30  40  M  W 


Cronk  Angle  (DEG) 


Nitric  Oxide  Concentration  Versus  Crank  Angle 


K1VA  Results  Tor  1500  ROM,  22  DBTDC  und  223  KP.  BMEP  Luud  (8  NM  Torque) 
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Cylinder  Gas  Pressure  Versus  Crank  Angle 


Cylinder  Gas  Temperature  Versus  Crank  Angle 


Mass  of  Fuel  Burnt  Venus  Crank  Angle 


Nitric  Oxide  Concentration  Versus  Crank  Angle 


KIVA  Results 


for  1500  RFM,  22  DBTDC  aad  335  KPa  BMEF  Load  (12  NM  Torqoe) 
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Maas  of  Burnt  Fuel  (G)  Goa  Pressure  (ATM) 


Crank  Angle  (DEG) 


Crank  Angle  (DEG) 


Cylinder  Gas  Pressure  Versus  Crank  Angle  Cylinder  Ga,  Temperature  Versus  Crank  Angle 


Crank  Angle  (DEG) 


Mass  of  Fuel  Burnt  Versus  Crank  Angle 


Nitric  Oxide  Concentration  Versus  Crank  Angle 


K1VA  Results  for  1500  RPM,  22  DBTDC  and  391  KPa  BMEP  Load  (14  NM  Torque) 
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Diesel  Combustion  Modeling  at  the  Engine  Research  Center 

C.  J.  Rutland  and  R.  D.  Reitz 
Engine  Research  Center 
Department  of  Mechanical  Engineering 
University  of  Wisconsin  -  Madison 
Madison,  WI 53706 

Modeling  diesel  combustion  presents  several  unique  challenges.  In  addition  to  turbulence  effects 
that  must  be  considered,  the  combustion  occurs  in  several  distinct  phases:  ignition,  premixed,  and 
diffusion.  The  combustion  models  also  rely  heavily  on  accurate  spray  breakup  and  evaporation 
models  to  predict  the  gaseous  fuel  distributions  prior  to  combustion. 

At  the  ERC  two  major  combustion  modeling  approaches  are  being  pursed.  The  first  is  a 
characteristic  time  scale  based  model.  This  approach  has  been  under  development  for  some  time 
and  is  performing  very  well  for  a  variety  of  engine  conditions.  The  second  approach  is  based  on 
flamelet  concepts  and  attempts  to  model  the  physics  more  directly.  This  method  is  newer  and  still 
under  development,  but  is  performing  fairly  well  for  several  cases.  The  flamelet  approach  offers 
more  potential  as  turbulence  modeling  moves  towards  large  eddy  simulations  and  as  emissions 
models  are  coupled  into  combustion  models. 

Both  approaches  use  ignition  modeling  based  on  the  Shell  low  temperature  multi-step  kinetics 
model.  This  is  a  reduced  mechanism  consisting  of  eight  reactions  using  fuel,  oxygen  and  five 
generic,  intermediate  species.  The  model  was  designed  to  simulate  auto-ignition  and  has  been 
adjusted  to  work  well  for  diesel  ignition  over  a  wide  range  of  conditions. 

In  the  time  scale  based  combustion  model  the  reaction  rate  is  determined  by  the  product  of  an 
inverse  time  scale  and  the  difference  between  the  local  species  mass  fractions  and  their  local 
thermodynamic  equilibrium  values.  The  time  scale  is  a  linear  combination  of  a  global,  Arrhenius 
reaction  time  scale  and  a  turbulent  mixing  time  scale.  The  Arrhenius  time  scale  represents  the 
laminar  burning  and  the  mixing  time  scale  represents  the  turbulent,  diffusion  combustion.  An 
exponential  sliding  factor  is  used  to  smoothly  transition  between  the  time  scales. 

The  flamelet  model  uses  a  similar  global  Arrhenius  reaction  for  the  premixed  combustion. 
However,  for  the  turbulent  diffusion  combustion,  a  coherent  flamelet  model  is  used.  In  this 
approach,  the  reaction  rate  is  the  product  of  a  flame  area  per  unit  volume  and  a  laminar,  flamelet 
consumption  rate.  The  flame  area  is  determined  from  a  transport  equation  that  includes  turbulent 
mixing  and  terms  that  represent  flame  area  sources  and  sinks.  The  main  source  of  flame  area  is 
turbulent  stretching.  Sink  or  destruction  mechanisms  are  included  for  flame  sheet  collisions, 
reactant  depletion,  front  propagation,  and  wall  effects. 

Both  combustion  modeling  approaches  require  criteria  and  mechanisms  for  modeling  the 
transition  between  the  different  types  of  burning.  The  main  criteria  for  the  end  of  ignition  and  start 
of  premixed  burning  is  the  local  temperature.  In  the  flamelet  model  this  is  augmented  by  a  second 
criterion  based  on  heat  release  rates. 

The  transition  between  premixed  combustion  and  turbulent  diffusion  burning  is  more  difficult. 
In  the  time  scale  model  it  is  handled  by  monitoring  the  amount  of  local  combustion  products  to 
indicate  the  progress  of  the  reaction.  In  the  flamelet  model,  the  transition  uses  a  Damkohler  number 
criterion  and  occurs  when  the  mixing  time  is  slow  compared  to  the  reaction  time  scale.  Then, 
premixed  burning  can  continue  with  previously  vaporized  fuel,  but  as  the  remaining  fuel  vaporizes 
it  goes  into  the  turbulent  diffusion  model. 

Overall,  the  combustion  models  compare  very  well  with  engine  data.  The  flamelet  model  has 
shown  good  agreement  on  several  cases  to  date  and  is  able  to  accurately  capture  the  different  types 
of  combustion.  The  time  scale  model  has  shown  excellent  agreement  for  a  wide  range  of  injection 
timings  and  pressures,  and  for  many  difficult  to  model  split  injection  cases. 


Work  supported  by:  ARO  Contract  #  DAAL03-92-G-0122,  DOE/NASA  Lewis  Contract  #  NAG3- 
1087,  Caterpillar;  with  computational  resources  provided  by  TACOM,  Cray  Research,  and  SDSC. 
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Diesel  Combustion  Modeling 
at  the  Engine  Research  Center 


C.  J.  Rutland 
R.D.  Reitz 


Engine  Research  Center  University  of  Wisconsin  -  Madison 

Work  Supported  by:  ARO  Contract # DAAL03-92-G-0122, 

DOE/NASA  Lewis  Contract  #  NAG3-1087 
Caterpillar 

Computational  Support  from:  TACOM,  Cray  Research,  and  SDSC 


Diesel  Combustion 


kinetics  controlled 


mixing  controlled 


Spray 

-  Atomization 

-  Breakup 

-  Vaporization 

Combustion 
•  Ignition 

-  Premixed  burning 

-  Diffusion  burning 
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Spray  Modeling 


•  Atomization  and  breakup  models 

-  Use  unstable  wave  growth  analysis 

-  Provides  breakup  time  scale  and  droplet  size 

•  Drag  modified  to  include  droplet  distortion 

•  Breakup  parameters  adjusted  for  wall 
interactions 

•  Multi-component  vaporization  completed  and 
currently  being  tested 


Combustion  Modeling  Approaches 


•  Characteristic  Time  Scale  Approach 

-  Very  well  developed 

-  Performing  very  well  over  a  wide  range 
of  engine  conditions 

•  Flamelet  Approach 
-Newer  model 

-  Good  performance  on  several  cases 

-  Offers  more  potential: 

-  Better  fit  with  large  eddy  simulations 

-  Better  coupling  to  emissions  models 
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Ignition  Modeling 


•  Based  on  Shell  Model 

-  Low  temperature,  multi-step  kinetics 

-  Uses  Fuel,  Oxygen  and  5  generic  species 

-  Has  8  reactions  representing  initiation, 
branching,  propagation  and  termination 

•  Modified  for  Diesel  Ignition 

-  Validated  using  bomb  data  for  diesel  fuel 

-  Applied  independently  in  each  grid  cell 


Characteristic  Time  Scale 
Combustion  Modeling 


•  Uses  Local  and  Equilibrium  Concentrations 

_  Local  concentration  -  Equilibrium  concentration 
X  laminar  +  /  X  turbulent 

•  Uses  Characteristic  Laminar  and  Turbulent 
Time  Scales 

X  laminar  =  A  exp(E/RT) 

X  turbulent  =  C2  k  /  e  (f  -  delay  coefficient ) 
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Flamelet  Modeling 


•  Premixed  Combustion 

-  Uses  same  laminar  global  reaction  model  as 
the  characteristic  time  scale  approach 

-  Premixed  burning  can  continue  in  same  cell 
in  which  diffusion  combustion  has  begun 

•  Diffusion  Combustion 

-  Based  on  coherent  flamelet  approach 

-  Reaction  rate  uses  local  chemistry  and 
flame  area  per  unit  volume 

Rate  =  p  VD  H 


Diffusion  Combustion  Model 


•  Laminar  Mass  Burning  Rate:  p  VD 

-  From  ID  analysis  of  diffusion  flame 
•  Uses  modified  local  concentrations 

•  Flame  Area:  I 

-  Obtained  from  turbulent  transport  equation 

|(pi)+V.  ((pi)u) =V.^vH+  apesi-Pf(i2) 

-  Source  for  I  from  turbulent  stretching 

-  Destruction  of  I  from: 

reactant  depletion  wall  effects 

flame  interaction  over  stretching 
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Diesel  Flamelet  Modeling 


•  Vaporizing  Fuel 

-  initially,  the  vaporizing  fuel  goes  to 
premixed  burning 

-  After  diffusion  flame  is  established  in  a  cell, 
vaporizing  fuel  goes  to  diffusion  burning 

•  Implementation 

-  Transport  equation  for  vaporized  premixed 
fuel  mass  fraction 

-  Transport  equation  for  inerts  to  determine 
actual  flame  stoichiometry 


Physical  Interpretation 


• 

m 

Ignition 

Premixed  Burn 

Diffusion  Burn 

-  Kinetic  controlled 

-  Shell  model 

•  Kinetic  controlled 

-  Global  Arrhenius 
model 

•  Mixing  controlled 

•  Flamelet  model 

-  Premixed  burn 
continues  with 
old  vaporized  fuel 
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Ignition  to  Premixed  Transition 


•  Local  Temperature  Criteria 

-  Switch  to  high  temperature  kinetics 
(premixed  combustion)  when  T  >  1000  K 

•  Used  by  both  modeling  approaches 

•  Heat  Release  Rate  Criteria  (fiamelet  approach) 

-  Use  linear  combination  of  ignition  and 

premixed  heat  release  rates  during  transition 

■  :•  /■ 

-  Helps  prevents  instabilities  during  transition 


Premixed  to  Diffusion  Transition 


•  Characteristic  Time  Scale  Approach 

-  Combine  laminar  and  turbulent  time  scales  with 
exponential  sliding  factor,  / 

-  Transition  factor  increases  as  reaction  progresses 

-  Reaction  progress  monitored  by  products 

•  Fiamelet  Approach 

-  Mixing  controlled:  diffusion  flame  initiated  when 
local  Damkohler  number  becomes  large 

Da  Ae*v{-EJRT) 

_ efk 

-  Initial  flame  area  equal  to  local  drop  surface  area 
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Engine  Simulations 


Caterpillar  Engine 

single  cylinder  3406 


-  Compression  ratio  14.5 

-  Displacement  2.44  liters 

-  Injector  6  holes,  up  to  90  MPa 

-  Speed  and  equiv  ratio  1600  rpm,  0.48 


Pressure  Traces 


Flamelet 


-50  0  50 

Crank  Angle  Degrees  ATDC 


Cat  Engine  1600  rpm 
SOI  -13  <t>  =  0.48 


Computed 

Measured 
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Heat  Release  Rate  (KJ/s) 


Heat  Release 


Characteristic  Time  Scale 


Flamelet 


Cat  Engine  1600  rpm 
SOI  -13  O  =  0.48 
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Heat  Release  Contours:  Flamelet 


•4  Crank  Angle  Degree  ATDC 


Heat  Release  Contours:  Flamelet 


1  Crank  Angle  Degree  ATDC 
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Detailed  Results  from  the 
Characteristic  Time  Scale  Model 


•  Effects  of  Injection  Timing 

•  Effects  of  Injection  Pressure 

•  Results  of  Emissions  Modeling 

•  Split  Injection  Results 

•  Alternative  Turbulence  Modeling 


Effects  of  Injection  Timing 


Computed 

Measured 


Crank  angle  (degrees) 


Effects  of  Injection  Pressure 


Injection  Pressure  and  Emissions 


NOx  (g/bhp-hr) 
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Emissions  Modeling 


Temperature 

H=  2370  K 
L=  958  K 


NOx-Mass  Fraction 

H=0.60  g/kg 
L=0.06  g/kg 


Soot-Mass  Fraction 

H=0.80  g/kg 
L=0.08  g/kg 

Cat  Engine  7°ATDC 


Split  Injection  Results 


3  degree  dwell 

- Computed  | 

-  Measured 

i 

1 
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Alternative  Turbulence  Model 

Results  for  RNG  modification  of  e  equation 


Summary 

•  Two  Approaches  to  Combustion  Modeling 

-  Characteristic  time  scale  and  flamelet  methods 

-  Both  methods  perform  well 

•  Modeling  Provides  Important  Information 

-  Effects  of  injection  pressure  and  timing 

-  Detailed  information  about  flow,  heat  release 
and  emissions 

•  Future  Directions 

-  Improve  burning  transitions 

-  Improve  spray  models 

-  Improve  turbulence  and  spray  mixing  models 
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